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ABSTRACT 

Liquid  crystals,  a  material  that  exhibits  brilliant  changes  in 
color  over  narrow  temperature  bands,  have  been  successfully  used  to 
study  the  temperature  field  produced  by  spherical  and  cylindrical 
radio  frequency  surgical  probes.  An  agar-water  mixture  was  used  to 
simulate  tissue.  Experimental  temperature  data  were  obtained  by 
photographing  the  liquid  crystal  material  through  the  clear  agar. 
Comparison  of  experimental  data  with  an  existing  theoretical  solution 
was  excellent,  agreeing  within  8%   for  the  spherical  probe.  However, 
the  cylindrical  data  and  cylindrical  theory  agreed  only  to  within  30%, 
the  large  discrepancy  being  attributed  to  edge  effects. 
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NOMENCLATURE 


Symbol 

Description 

Typical  Units 

A 

Surface  area 

cm2 

o 

A 

Wavelength 

Angstrom 

b 

Number  of  flow  tubes 

C 

Volumetric  heat  capacity 

cal/cm3  °( 

°c 

Temperature  in  degrees  Centigrade 

D 

Diameter 

cm 

E 

Electric  potential 

volts 

f 

Frequency 

hertz 

G 

Nondimensional  power 

G' 

Corrected  nondimensional  power 

I 

Current  density 

amps/cm2 

I 

Radio  frequency  current 

amps 

i 

Coordinate  of  nodal  point  along 
distance  axis  in  finite 
difference  formulation 

k 

Thermal  conductivity 

W/m  °K 

L 

Length 

cm 

M 

Stability  criterion  for  finite 

difference  formulation 

m       Slope  of  nondimensional  resistivity 
curve 

n       Coordinate  of  nodal  point  along 
time  axis  in  finite  difference 
formulation 

R       Nondimensional  radius  to  an  isotherm 

R-jq      Ten  ohm-ten  watt  resistor 


r 

0 

Probe  radius 

r 

Radius 

rt 

Isotherm  radius 

T 

Temperature 

t 

Time 

V 

Voltage 

a 

Thermal  diffusivity 

cm 

cm 

cm 

°C 

seconds 

volts 

cm2/sec 

Y      Slope  of  dimensional  resistivity 
curve 

n      Ratio  of  speed  of  light  in  air 
to  that  in  water 

e  Nondimensional  temperature 

A  Wavelength  Angstrom 

p  Resistivity  ohm- cm 

p  Nondimensional  resistivity 

t  Nondimensional  time 

<f>  General  potential 

to  Uncertainty 

v  Del  operator 

SUBSCRIPTS 

o  Initial  conditions 

10  At  the  ten  ohm  resistor 

cct  In  the  circuit 

i  Of  the  ith  component 

p  From  the  photographs 

cyl  With  reference  to  the  cylindrical  probe 

sphere  With  reference  to  the  spherical  probe 
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NONDIMENSIONAL  QUANTITIES 
T  -  T 


D 

T 

0 

T 

= 

at 

r2 

0 

I2p    r2 

G 

0 

= 

0      0 

A^  k  T 

r 

0 


Temperature 


Time 


Power 


Distance 


p   -  -~-    =  1  +  me        Resistivity 

Po 
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I.   INTRODUCTION 

In  neurosurgery,  radio  frequency  current  has  been  used  to  create 
lesions,  thus  alleviating  conditions  such  as  acute  pain  and  Parkinsonism. 
A  current  emitting  electrode,  similar  to  a  hypodermic  needle,  is  inserted 
into  the  brain  through  the  skull.  An  inactive  electrode  is  placed  else- 
where on  the  body  in  an  area  of  high  conductivity.  Radio  frequency 
current,  at  a  low  power  level,  is  passed  between  the  electrodes  and 
electrical  energy  is  transformed  to  heat  in  the  tissue  with  the  result- 
ant rise  in  temperature  destroying  the  tissue. 

To  accurately  use  the  radio  frequency  technique,  a  surgeon  must 
know  the  exact  location  of  the  electrode  (probe)  tip  in  the  brain,  the 
volume  of  tissue  to  be  destroyed,  the  thermal  and  electrical  properties 
of  the  tissue  [1],  and  the  temperature  near  the  probe  tip.  Thermocouples 
and  thermistors  attached  to  the  .probe  have  been  used  to  monitor  probe  tip 
temperatures  [2,  3,  4].  However,  these  temperature  sensors  indicated  the 
temperature  of  the  probe,  not  the  tissue  temperature  itself.  Addition- 
ally, the  sensors  were  nearly  as  large  as  the  probe  tip  and  may  have 
distorted  the  thermal  field  by  acting  as  thermal  sinks.  Without  accurate 
knowledge  of  tissue  temperature,  tissue  boiling  occurred  with  attendant 
gas  formation  and  tissue  carbonization  [4].  These  effects  resulted  in 
irregular  shaped  lesions  and  unpredictable  damage  along  the  entrance 
path  of  the  probe  into  the  tissue. 

The  objectives  of  this  thesis  were  twofold.  First,  to  experimentally 
determine  the  thermal  field  produced  by  a  radio  frequency  probe  embedded 
in  a  resistive  agar  medium.  Liquid  crystals,  a  material  that  indicates 
temperature  through  a  color  change,  were  used  as  the  temperature  sensors. 
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The  second  objective  was  to  compare  experimental  data  with  results  pre- 
dicted using  a  theoretical  heat  transfer  model  of  the  radio  frequency 
probe  [5],  The  general  theoretical  model  included  the  effects  of  tissue 
metabolism  and  blood  profusion.  These  two  effects  have  been  excluded 
from  the  present  investigation  since  the  thermal  field  was  generated  in 
an  inert  gel-like  substance. 

The  temperature  fields  produced  by  spherical  and  cylindrical  radio 
frequency  probes  were  successfully  studied  with  the  liquid  crystal 
material.  ,  The  crystals,  with  their  low  thermal  capacity  [6],  produced 
brilliant,  undistorted,  continuous  displays  of  several  selected  isotherms 
surrounding  the  probes.  Agreement  between  the  theoretically  predicted 
and  experimentally  measured  values  of  temperature  for  the  spherical 
probe  was  excellent  (8%).  Agreement  between  theory  and  experiment  for 
the  cylindrical  probe  was  not  good  (30%).  The  large  differences  observed 
in  the  cylindrical  case  are  attributed  to  edge  effects.  The  theory  assumed 
an  infinitely  long  probe  while  the  experiment  employed  probes  with  length 
to  diameter  ratios  of  4,  10,  and  20.  The  edge  effects  were  clearly  dis- 
played by  the  liquid  crystals. 
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II.  HISTORICAL  BACKGROUND 

Use  of  radio  frequency  current  to  create  physiologic  lesions  is  not 
new.  Carpenter  and  Whittier  [7]  describe  early  investigations  into  the 
uses  of  radio  frequency  current  as  a  surgical  tool.  Aronow  [1]  was  the 
first  investigator  of  recent  time  to  study  the  problem  of  radio  frequency 
current  in  tissue.  Lesion  formation  using  radio  frequency  current  is 
accomplished  by  inserting  into  the  tissue  a  small  hypodermic  needle  which 
is  electrically  insulated  except  at  the  tip.  This  electrode  acts  as  the 
active,  current  emitting  electrode.  A  larger,  inactive  electrode  is 
placed  elsewhere  on  the  body.  The  path  of  current  conduction,  i.e.  the 
body,  acts  as  a  leaky  dielectric  by  dissipating  electrical  energy  in  the 
form  of  heat  [1].  This  Joulean  heating  effect  results  in  a  local  rise  in 
tissue  temperature.  When  tissue  temperature  is  increased  to  55°C,  tissue 
destruction  occurs.  Accurate  prediction  of  when  tissue  in  the  probe  tip 
region  has  reached  this  temperature  has  been  a  major  problem  with  the 
radio  frequency  probe. 

Excessive  tip  temperature  causes  gas  formation  and  tissue  carboniza- 
tion which  results  in  poor  prediction  of  lesion  size  and  additional 
damage  along  the  path  of  electrode  insertion  [4].  Studies  with  thermo- 
couples and  thermistors  have  been  conducted  to  correlate  electrode 
temperature  with  lesion  size  [2,  3].  Since  heating  takes  place  in  the 
tissue  and  not  in  the  electrode,  the  sensors  inside  the  probe  tip  indi- 
cate only  that  temperature  rise  caused  by  conduction  from  the  tissue  back 
into  the  electrode.  This  temperature  rise  may  be  less  than  the  maximum 
tissue  temperature.  Therefore,  reliable  prevention  of  excessive  tip 
temperature  has  been  difficult.  The  same  type  of  tissue  temperature 

measurement  problem  occurs  with  the  direct  current  probe,  [_5j. 
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The  direct  current  probe  has  a  more  serious  liability  than  the 
tissue  temperature  measurement  problem.  Direct  current  has  preferred 
paths  in  the  body  which  results  in  irregular  lesions.  Death  may  occur 
if  the  direct  current  finds  a  path  to  the  heart.  Preferred  paths  in 
the  body  do  not  appear  to  be  a  problem  with  the  radio  frequency  or  the 
high  resistance  probe. 

The  high  resistance  probe  consists  of  a  high  resistance  iron-nickel 
alloy  wrapped  around  a  nonconductive  rod.  Electrical  current  is  passed 
through  the  wire  causing  the  probe  to  be  heated.  Groff's  study  of  the 
high  resistance  probe  [8]  demonstrated  the  feasibility  of  producing 
controlled,  predictable  heated  regions  in  agar.  In  opposition  to  the 
radio  frequency  probe,  where  heating  occurs  in  the  tissue,  the  high 
resistance  probe  itself  generates  heat.  Thermal  energy  is  then  trans- 
ferred to  the  tissue  by  conduction.  Prediction  of  the  magnitude  of  the 
temperature  field  produced  by  a  high  resistance  probe  thus  requires 
knowledge  of. only  the  thermal  properties  of  the  medium.  Another  advan- 
tage of  the  high  resistance  probe  is  the  fact  that  a  surgeon  may  check 
for  proper  probe  tip  location  before  forming  the  lesion. 

Probe  tip  location  may  be  verified  by  raising  the  tissue  temperature 
to  40-49°C  using  low  current  levels.  Increasing  temperature  into  this 
range  blocks  nerve  functioning  but  does  not  cause  permanent  damage.  This 
technique,  equally  applicable  to  the  high  resistance  and  the  radio  fre- 
quency probe,  is  especially  helpful  in  treating  Parkinsonism.  As  the 
temperature  increases,  shaking  diminishes  and  stops  when  the  probe  tip 
has  been  properly  located.  Lesion  formation  is  then  accomplished  by 
increasing  the  local  tissue  temperature  to  55°C,  being  careful  to  prevent 
excessive  tip  temperatures  that  result  in  tissue  boiling. 
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While  the  high  resistance  probe  requires  knowledge  of  only  the 
thermal  properties  of  tissue,  the  radio  frequency  probe  requires 
knowledge  of  both  thermal  and  electrical  properties.  Geddes  and 
Baker  [9]  discuss  difficulties  in  obtaining  electrical  properties  of 
tissue.  Since  tissues  are  composed  of  cells,  different  properties 
may  be  exhibited  depending  on  the  cell  orientation.  Other  factors 
influencing  the  electrical  property  values  include  age,  history,  and 
environment  of  the  tissue.  Tissue  properties  measured  IN  VITRO 
differ  from  those  measured  IN  VIVO.  Thermal  properties,  on  the  other 
hand,  are  reliably  predicted  based  on  the  tissue  water  content  [10], 
The  prediction  of  lesion  size  based  on  the  magnitude  and  duration  of 
radio  frequency  current  and  tissue  electrical  properties  is,  there- 
fore, more  complex  when  using  the  radio  frequency  probe  than  when 
using  the  high  resistance  probe. 
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III.  THEORETICAL  DEVELOPMENT 

A.  SPHERICAL  PROBE 

Consider  a  spherical  probe,  of  radius  r  ,  which  continuously  emits 
radio  frequency  current,  I,  into  a  medium  that  is  homogeneous  and  infinite 
in  extent.  Assume  that  the  medium  possesses  electrical  resistivity,  p, 
thermal  conductivity,  k,  and  volumetric  heat  capacity,  C.  Gengler  [5] 
has  derived  the  basic  equation  which  accounts  for  the  thermal  effects 
produced  in  the  medium  due  to  the  Joulean  heating  effects  of  the  radio 
frequency  current: 


r  hi 

iLkr2  3rJ 


1   3lkr2  3rJ     .  I2P     r  3T  ,,  N 


In  the  present  investigation,  the  experimental  test  medium,  an  agar- 
water  mixture,  was  found  to  exhibit  an  electrical  resistivity  that 
varied  linearly  with  temperature: 


P  =  P0  +  t(T  -  T  ) 


(2) 


where  the  subscript  indicates  an  initial,  reference  condition. 
Substitution  of  Equation  (2)  into  Equation  (1)  yields: 

T        all 

1       alkr2  3rJ         I2P  I2y   (T  -  To)  3T 

T7         ar  16Tr2r4    +      Tefir*  "     L  3t 


(3) 


Assuming  that  k  and  C  are  constant,  Equation  (3)  may  be  normalized  by 
introducing  the  following  non-dimensional  quantities: 
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R 

= 

r 
r 

0 

Distance 

T 

= 

kt 
Cr  2  " 

0 

at 

r  2 

0 

Time 

0 

= 

T  -  T 

0 

T 

0 

Temperature 

G 

0 

= 

I2P  r2 

0  0 

Power 

k  T 

0 

P 
m 

= 

p    =  1 

po 

Po 

+  me 

Resistivity 

where 

7  = 

I 

I 

"  A 

the  current 

4irr2 

0 

Substitution  into  Equation   (3)  yields 


6 
^R~  '  R^ 


+  ^  (l  +  me)   -    §i 


(4) 


The  boundary  and  initial  conditions  are: 

a)  at  the  probe  surface,  there  is  no  heat  conduction  back  into 
the  probe, 

|ie.|      =  0;  (5) 

\9R/  R  =  1 

b)  the  temperature  field  is  undistorted  at  large  distances 
from  the  probe, 

R  -+  »,  e  +  0;  (6) 
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c)  the  medium  is  initially  at  uniform  temperature, 

t  =  0,  e  =  0.  (7) 

Equations  (4),  (5),  (6),  and  (7)  may  be  cast  in  explicit  finite  dif- 
ference form  as  [5] 


aR< 


1   -  4- 


f  K-i +     +  tK+i  " 2 


TW 


(4a) 


mG  e. 

0     1 

(iAR)4 


>r  -  j 

Ax 


n  n 

'i*_l   =  V+l 


(5a) 


n  n 

>N-1   "  eN+l 


(6a) 


ei  =  0.0 


(7a) 


where  i  and  n  define  the  coordinates  of  a  nodal  point  in  an  R,  x 
network  such  that  the  nondimensional  temperature  at  time  t  =  nAx 
and  radial  location  R  =  iAR  is  defined  as  e"  ;  i*  is  at  R  =  1; 
i  =  i*,  i*  +  1,...,N-1,N.  Solving  Equation  (4a)  for  e.   yields: 
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n+1 


At 
AR^ 


tKH1+TK+1 


,At 


AxmG     \     AxG 
o    1 .  o 


AR7    +    JurJ"  +  (iAR)4 


Define  M  =  At/aR2,  then 


n+1 


=  M 


1  - 


'?-! 


1   + 


>?+1 


ev  I   -  2M  +  - 


MmAR2Go]   AxG0 
(iART^j  (iAR)4 


Since  an  explicit  finite  difference  scheme  has  been  employed,  a  stability 
criterion  must  be  stated  to  limit  the  size  of  the  time  step  based  on  the 
value  of  incremental  distance  (aR).  Discussions  of  stability  may  be  found 
in  References  11  and  12.  The  stability  criterion  for  this  problem  becomes 


MmAR2G 


1  -  2M  -» 


TTarJ 


T~ 


>  0 


Solving  for  M  and  assuming  the  worse  case, 


M  = 


1 


(9) 


2  -  n;AR2(?p 

c       (iAR)^ 


Since  M  was  originally  defined  as  At/aR2,  choosing  a  value  of  aR2  fixes 
the  value  of  At. 

Computer  program  I  (Appendix  E)  was  used  to  solve  Equation  (8)  subject 
to  boundary,  initial,  and  stability  conditions  (Equations  5a,  6a,  7a,  and 
9).  Reference  13  states  that  the  error  in  an  explicit  difference  approxi- 
mation is  on  the  order  of  aR2.  By  varying  the  size  of  the  aR  increments 
and  plotting  temperature  versus  aR2,  it  was  possible  to  extrapolate  the 
exact  value  of  temperature.  Figure  1  represents  this  technique,  known  as 
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Richardson's  technique,  for  the  spherical  probe.  It  can  be  seen 
that  using  a  value  of  aR  =  0.1  yielded  theoretical  results  within 
2%   of  the  exact  value. 

B.   CYLINDRICAL  PROBE 

Paralleling  the  development  of  the  basic  equations  for  the  spheri- 
cal probe,  the  equation  for  the  cylindrical  probe  in  nondimensional 
form  is: 
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Substituting  G  =  G  (1  +  me)  into  Equation  (10)  and  transforming  to 
finite  difference  form  yields: 
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The  appropriate  boundary  conditions  are 
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where  i*  corresponds  to  R  ■  1.0 

i  =  i*,  1*  +  1,...,  N-1,  N 
The  stability  criterion  for  the  cylindrical  case  is: 

M  B      1 

mG 

2  -  -j2— 

where  M  is  defined  as  in  the  spherical  case  as  At/aR2. 

Computer  program  II  (Appendix  E)  was  used  to  generate  the  numerical 

solution  to  this  problem. 
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IV.  EXPERIMENT 

A.   EQUIPMENT 

A  radio  frequency  generator,  Model  RFG-2AV,  manufactured  by 
Radionics  Incorporated  of  Burlington,  Massachusetts,  was  used  in  the 
present  investigation  to  generate  radio  frequency  current.  Output 
from  the  generator  was  a  continuous,  unmodulated  sine  wave  at  500 
kilohertz  with  a  maximum  power  of  18  watts.  Although  the  generator 
included  a  built-in  radio  frequency  voltmeter  and  radio  frequency 
milliammeter,  two  independent  RMS  meters  were  utilized  in  data  collec- 
tion to  allow  the  use  of  more  sensitive  scales  where  required.  Figure  2 
displays  the  generator,  voltmeters,  and  one  of  the  experimental  probe 
test  cells  used  in  the  experiments.  Figure  3  is  a  schematic  of  the 
circuitry  used. 

Both  spherical  and  cylindrical  radio  frequency  probe  geometries 
were  studied  experimentally.  An  agar-water  mixture  (a  gel -like 
substance)  was  used  to  simulate  tissue.  The  mixture  was  contained  in 
the  test  cells  shown  in  Figure  4.  In  both  cells,  at  least  25  probe 
radii  separated  the  probe  from  the  nearest  solid  boundary.  Thus  the 
cells  could  be  considered  as  surfaces  at  infinity  in  relation  to  the 
probes.  Plexiglass  was  used  for  the  probe-agar  cells  due  to  its  trans- 
parence and  because  it  is  a  dielectric  at  the  low  power  levels  used 
during  experimentation.  Boundary  conditions  were  changed  from  insulated 
surfaces  to  conductive  surfaces  by  lining  the  inside  surfaces  of  the 
cells  with  conductive  foil.  The  current  path  was  then  from  the  probe, 
through  the  agar  medium,  to  the  foil  with  electrical  energy  dissipated 
as  heat  in  the  anar. 
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PROBE   CELL 


RADIO  FREQUENCY 
GENERATOR 


Voltmeter  V-,  measured  the  voltage  drop  across  the 
ten  ohm  resistor  R1Q  .  Circuit  current  was  deter- 
mined from  the  two  quantities  using  Kirchoff  s  Law. 


Circuitry  for  radio  frequency  probe  experiments 


Figure  3 
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Local  temperature  increases  in  the  agar  were  measured  by  use  of 
liquid  crystals.  Appendix  B  presents  a  complete  discussion  of  liquid 
crystals.  Briefly,  liquid  crystals  change  color  at  specific  temperatures, 
the  particular  temperature  depending  on  the  crystal  in  use.  The  color 
change  thus  allows  a  continuous  display  of  isotherm  fronts  for  the  map- 
ping of  thermal  regions. 

In  the  experiments,  liquid  crystals  manufactured  by  National  Cash 
Register  and  Hoffman-LaRoche  were  used.  The  Hoffman-LaRoche  liquid 
crystal,  ROCHROME,  was  provided  in  the  form  of  a  tape  which  could  be 
applied  directly  to  the  object  under  study.  The  National  Cash  Register 
liquid  crystals  were  obtained  in  liquid  form.  The  crystals  are  actually 
encapsulated  in  spheres  5  to  30  microns  in  diameter,  thus  preventing 
contamination  of  the  crystal.  Due  to  this  encapsulation  procedure,  it 
was  possible  to  mix  different  liquid  crystals  together  without  changing 
the  characteristics  of  any  single  crystal.  Several  distinct  sets  of 
isotherms  were  viewed  since  the  liquid  crystals  mixed  did  not  have  tem- 
perature ranges  that  overlapped. 

In  the  experiments,  three  liquid  crystals  were  mixed  together  and 
were  applied  to  a  transparent  0.003  inch  mylar  substrate.  The  liquid 
crystals  were  then  covered  with  black  paint  and  polyurathane  varnish. 
The  black  background  was  required  to  absorb  the  light  transmitted 
through  the  liquid  crystals;  the  polyurathane  varnish  sealed  the  black 
paint.  Sealing  the  system  in  this  manner  allowed  repeated  use  of  the 
liquid  crystal  system.  Figure  5  displays  the  liquid  crystal  sheet 
encased  in  the  holders  used  in  the  test  cells. 

Prior  to  experimental  trials,  the  crystals  were  calibrated  using 
the  Rosemount  Constant  Temperature  Bath  following  the  procedure  of 
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Petrovic  [14].  The  following  liquid  crystals  were  calibrated  and 
subsequently  used  in  experimental  trials: 

Table  1 
Manufacturer  Designation    Red      Green    Blue 

National  Cash  Register  R-27         26.7°C    27.9°C    28.9°C 
Hoffman-LaRoche        Rochrome      29.9°C    31.5°C    34.2°C 
National  Cash  Register  R-37         35.7°C    37.0°C    38.5°C 
National  Cash  Register  R-53         50.5°C    52.TC    53.3°C 
Data  were  collected  at  specified  times  after  application  of  power  by 
photographing  the  isotherms  displayed  on  the  liquid  crystal  sheet. 
Distances  from  the  probe  center  to  an  isotherm  (color  band)  under  consid- 
eration were  easily  measured  from  the  photographs. 

B.   EXPERIMENTAL  PROCEDURE 

The  agar  solution  was  prepared  by  adding  3.0  grams  of  DIFCO 
"BACTO-AGAR"  powder  to  997  grams  of  boiling  distilled  water.  When  the 
powdered  agar  had  thoroughly  dissolved,  a  100  ml.  sample  was  taken  to 
fill  the  resistivity  test  cell  which  was  used  to  determine  the  variation 
of  resistivity  with  temperature  (see  Appendix  A).  Prior  to  jelling,  the 
remaining  agar  solution  was  poured  into  one  of  the  probe  test  cells.  The 
liquid  crystal  sheet  was  positioned  and  the  probe  located  such  that  a  probe 
diameter  was  in  the  plane  of  the  liquid  crystals.  The  entire  apparatus  was 
then  cooled  for  approximately  six  hours  to  achieve  a  uniform  temperature  in 
the  agar.  At  that  time  the  initial  agar  temperature  was  recorded,  the 
generator  circuit  was  connected,  and  power  was  applied.  Photographs  of 
the  liquid  crystal  sheet  were  taken  through  the  agar  at  specified  times, 
thus  recording  the  color  change  in  the  liquid  crystals.  Distances  were 
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measured  directly  from  the  pictures  and  were  scaled  with  respect  to 
the  probe  radius  to  form  a  nondimensional  distance. 

All  quantities  were  nondimensional ized  to  allow  general  application 
of  the  data.  The  following  nondimensional  quantities  were  used  (5): 
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r.  =   distance  from  the  probe  center  to  the  outer  edge  of  the 
isotherm  (color  band)  under  consideration  (cm) 

r   =   probe  rad.ius  (cm) 

a   =   thermal  diffusivity  (cm2/sec) 

t   =   time  (sec) 

T   =   isotherm  temperature  (°C) 

T   =   initial  agar  temperature  (°C) 

I   =   radio  frequency  current  (amps) 

p^  =   initial  agar  resistivity  (ohm-cm) 
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Using  a  constant  temperature-power  ratio  (e/G  ),  curves  were 
plotted  for  distance  (R)  versus  time  (t)  based  on  the  solution  of 
the  finite  difference  equations  for  the  probe  geometry  under  investv 
gation  (Equations  8  or  11).  The  curves  were  then  compared  to  the 
experimental  data  points  of  nondimensional  isotherm  location  (R)  at 
various  times  (x).  An  uncertainty  analysis  was  conducted  and  is 
presented  in  Appendix  D. 
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y.   EXPERIMENTAL  RESULTS 

A.  SPHERICAL  PROBE 

Figure  6  is  a  series  of  photographs  of  the  growth  of  the  radial 
temperature  field  around  the  spherical  probe  at  nondimensional  times 
of  1.5,  3.1,  6.2,  and  12.5  (corresponding  to  times  of  60,  120,  240, 
and  480  seconds  for  a  spherical  probe  of  radius  0.24  cm.).  In 
Figure  6c,  the  isotherms  used  in  data  collection  are  labeled  for 
reference.  The  temperature-color  relationships  used  in  data 
collection  were: 

Liquid  Crystal      Temperature        Color 

R-27  27.9°C  Green 

Rochrome  29.9°C  Red 

R-37  37.0°C  Green 

The  R-53  appeared  on  the  photograph  at  large  times,  but  was  not 
actually  used  in  data  reduction  and  comparison  with  theory.  Values  of 
nondimensional  radius  (R)  and  time  (t)  were  plotted  for  comparison  with 
the  theoretical  curves  based  on  temperature-power  ratios  (e/G  ). 
Figures  7  and  8  are  representative  of  the  comparisons.  The  dotted 
lines  are  the  uncertainty  limits,  which  are  discussed  in  Appendix  D. 
Experimental  results  were  within  8%  of  the  theoretical  values.  Thus,  the 
spherical  theory,  taking  into  account  variable  resistivity,  was  verified 
by  the  experimental  procedure  of  using  liquid  crystals  to  determine 
thermal  changes. 
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Figure  6a 
Tau  =  1.5   t  =  60  sec, 


Figure  6b 
Tau  =3.1   t  =  120  sec, 


Figure  6c 
Tau  =6.2   t  =  240  sec. 


Figure  6d 
Tau  =  12.5  t  =  480  sec, 


Figure  6.  Spherical  transient  temperature  field  growth 
with  probe  diameter  of  0.241  cm. 


33 


4.0 


3.0 
r, 


R=: 


l 


Figure  7 
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Figure  8 
Nondimensional  isotherm  radii 


vs. 


nondimensional  time  for  a  spherical 
probe  of  diameter  0.482  cm 
using  three  liquid  crystals. 
Data  set  J 
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B.   CYLINDRICAL  PROBE 

Figures  9,  10,  and  11  display  the  experimentally  determined 
transient  temperature  field  at  nondimensional  times  of  4.2,  6.25, 
9.3,  and  12.5  (corresponding  to  times  of  40,  60,  90,  and  120  seconds 
for  a  cylindrical  probe  of  radius  0.12  cm.).  The  length  to  diameter 
ratios  used  were  4,  10,  and  20  in  Figures  9,  10,  and  11,  respectively, 
Values  of  distance  (R)  were  plotted  as  a  function  of  time  (t)  for 
constant  temperature-power  ratios  (e/G  ),  as  seen  in  Figures  12,  13, 
and  14.  The  solid  lines  represent  the  theoretical  curves. 
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Figure  9a 
Tau  =  4.2   t  =  40  sec. 


Figure  9b 
Tau  =  6.2   t  =  60  sec, 


Figure  9c 
Tau  =9.3   t  =  90  sec. 


Figure  9d 
Tau  =  12.5  t  =  120  sec, 


Figure  9.  Cylindrical  transient  temperature  field 

growth,  L/D  =  4,  probe  diameter  =  0.241  cm. 
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Figure  10a 
Tau  =  4.2   t  =  40  sec. 


Figure  10b 
Tau  =  6.2   t  =  60  sec. 


Figure  10c 
Tau  =9.3   t  =  90  sec. 


Figure  lOd 
Tau  =  12.5  t  =  120  sec. 


Figure  10.  Cylindrical  transient  temperature  field 

growth,  L/D  =  10,  probe  diameter  =  0.241  cm. 
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Figure  11a 

4.2   t  =  40  sec, 


Figure  lib 
Tau  =  6.2   t  =  60  sec, 


Figure  lie 
Tau  =9.3   t  =  90  sec. 


Figure  lid 
Tau  =  12.5  t  =  120  sec, 


Figure  11.  Cylindrical  transient  temperature  field 

growth,  L/D  =  20,  probe  diameter  =  0.241  cm. 
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-  Nondimensional  isotherm  radii 
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Figure  13 

Nondimensional  isotherm  radii 

vs. 

nondimensional  time  for  a  cylindrical 

probe  with  L/D  =  10  and  diameter  =  0.241  cm. 
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Figure  14 
Nondimensional  isotherm  radii 
vs. 
nondimensional  time  for  a  cylindrical 
probe  with  L/D  =  20  and  diameter  =  0.241  cm. 
Data  set  K 
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VI.  DISCUSSION  OF  RESULTS 

A  sertous  difficulty  in  modeling  the  radio  frequency  probe  in 
tissue  has  previously  been  an  inability  to  accurately  map  the  thermal 
area  surrounding  the  probe.  Thermocouples  and  thermistors  have  either 
been  placed  inside  the  probe  tip,  where  they  actually  measured  the 
probe  temperature,  or  placed  in  the  medium,  where  they  acted  as  thermal 
sinks.  The  use  of  liquid  crystals,  as  demonstrated  in  this  thesis, 
now  allows  an  accurate,  continuous  display  of  specific  temperatures 
depending  on  the  crystal  in  use.  Liquid  crystals  provide  a  flexible 
tool  for  temperature  measurement  regardless  of  probe  geometry.  Figure  15 
shows  liquid  crystals  indicating  isotherms  which  were  caused  by  radio 
frequency  current  between  two  nonsymmetric  probes.  As  can  be  seen,  the 
liquid  crystals  react  to  provide  a  complete  map  of  the  thermal  region. 

In  the  spherical  geometry,  the  experimental  values  were  in  good 
agreement  with  the  theory  as  modified  for  variable  resistivity  (8%). 
The  cylindrical  data  and  theory  were  not  in  agreement  (30%).  As  a  result 
of  the  poor  agreement  between  theory  and  experiment  in  the  cylindrical 
geometry,  an  examination  of  the  flux  lines  emanating  from  the  probe 
was  conducted.  Appendix  C  discusses  the  flux  plots  obtained. 
Figure  16  is  a  flux  plot  showing  only  the  lines  of  constant  current. 
As  can  be  seen,  the  direction  of  the  current  lines  is  changed  in  the 
area  above  the  probe.  The  bending  of  the  current  lines  results  in  a 
higher  local  current  density  with  greater  heating  than  theoretically 
predicted. 
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Figure  15 

Liquid  crystals  map  a 

nonsymmetric  thermal  region 
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Figure  16.  Constant  current  lines  emanating  from 
a  cylindrical  probe  based  on  teledeltus  plot  of 
the  isopotentials  in  plane  perpendicular  to  probe 
axis 
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An  approximation  was  derived  to  account  for  the  increased  current 
density.  Due  to  symmetry  of  the  problem,  only  half  the  probe-cell 
system  was  modeled.  The  following  parameters  were  defined: 

Number  of  flow  tubes  (for  half  field)  =  b 

Current  flow  (for  half  field)  =  1/2 

Half  the  probe  surface  area  -  vr  I 

Angle  subtended  by  the  ith  flow  tube  =  <j>. 

Surface  area  of  the  ith  flow  tube  at  the 
probe  surface  =  <j>  .r  L 

1  I  0 

Current  flow  per  flow  tube  =  I/2b  =  constant 
Current  flux  in  the  ith  flow  tube  =  I/2b*.r  L 

1  0 

Recall  that  I  was  previously  defined  as 

I   =  -4-  ,  or  for  this  situation,  the  local  current  density  is 

i,  -,  ' 


i    2b<j).r  L 

1  o 

I.2p  r  2 
The  local  power  parameter  is  G.  =  — ■  T°  °     or  after  substitution 


for  1.2, 


o 


IT2 


1    o   §*  o 

From  Figure  16,  $  =  14°,  b  =  8  which  yields  G.  =  2.58  G  . 
Applying  this  correction  to  the  nondimensional  temperature-power  term 
(e/G  )  shifts  the  theoretical  curve  so  that  experimental  data  more 
closely  approximate  the  theory  (Figure  17). 

Since  the  degree  of  current  density  is  influenced  by  probe  depth, 
the  correction  factor  must  be  determined  for  each  experiment. 

Examination  of  Figure  18  demonstrates  the  strong  end  effects  on  the 
finite  cylindrical  probe.  The  cylindrical  theory  (Equation  11)  was 
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Figure  17 
Comparison  of  experimental  data  with 
original  theoretical  curve  and 
corrected  theoretical  curve 
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Figure  18.  Liquid  crystals  demonstrating  the 
pronounced  edge  effects  caused  by  a  cylindrical 
probe  (L'D  =  20) 
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based  on  an  infinite  cylindrical  probe  with  no  end  effects.  The  figures 
show  the  strong  thermal  fields  at  the  probe  ends  developing  before  the 
effects  at  the  probe  center  are  observed.  With  time,  the  conduction  of 
heat  in  the  agar  tended  to  reduce  the  thermal  gradient  between  the  probe 
ends  and  its  midpoint,  with  the  resultant  growth  as  shown  in  Figure  11. 
A  possible  explanation  for  the  early  observation  of  thermal  fields  at 
the  probe  ends  is  that  radio  frequency  current  density  is  higher  at 
sharp  corners  than  at  plane  surfaces.  The  higher  current  density  results 
in  greater  heating. 

By  using  a  probe  with  a  small  surface  area,  resulting  in  high 
current  density,  it  was  possible  to  observe  boiling  in  the  agar.  Since 
the  liquid  crystals  used  in  the  experiments  were  active  in  temperature 
ranges  below  the  boiling  temperature  of  the  agar,  no  data  were  taken  on 
boiling.  However,  Figure  19  is  a  photograph  taken  of  the  boiling 
phenomenon.   It  was  of  interest  to  observe  formation  of  a  bubble  as  the 
agar  was  heated  to  boiling.  When  rising,  the  bubble  traveled  until  it 
encountered  a  cool  region;  the  bubble  then  collapsed.  This  boiling  is 
analogous  to  that  occurring  in  actual  surgery.  The  reason  for  controlling 
the  probe  tip  temperature  is  readily  apparent  when  observing  the  deteri- 
oration of  the  jelled  agar  as  bubbles  rise  randomly  towards  the  surface; 
the  analogy  in  surgery  being  unpredictable  physiological  damage. 
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Figure  19.  Boiling  of  the  agar-water  mixture 
caused  by  excessive  probe  tip  temperature 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  use  of  liquid  crystals  for  mapping  regions  surrounding  heated 
surgical  probes  has  been  clearly  demonstrated.  When  accuracy  on  the 
order  of  1.0°C  is  sufficient ,  liquid  crystals  provide  a  continuous 
picture  of  thermal  changes  as  displayed  in  Figures  20,  21,  22,  and  23. 
Through  the  use  of  liquid  crystals,  the  spherical  theory  developed  by 
Gengler  [5]  has  been  verified.  On  the  other  hand,  Gengler's  cylindrical 
theory  has  been  shown  to  require  further  work  in  order  to  account  for  the 
finite  cylindrical  probe  with  its  strong  end  effects. 

It  is  most  probable  that  an  actual  radio  frequency  surgical  probe 
would  be  constructed  with  a  length  to  diameter  ratio  in  the  range  of 
two  to  five  [1,  2,  15,  16].  Clearly,  a  theoretical  treatment  for  an 
infinite  cylinder  would  break  down  at  these  small  length  to  diameter 
ratios.  As  a  demonstration  of  this  fact,  an  equivalent  surface  area 
for  a  spherical  probe  was  determined  based  on  the  surface  area  of  the 
cylindrical  probe  with  a  length  to  diameter  ratio  of  four.  Figure  24 
is  a  comparison  of  the  cylindrical  theory,  the  spherical  theory,  and 
the  experimental  data.  Relationships  between  the  cylindrical  and  the 
spherical  parameters  have  been  derived  based  on  the  fact  that  the  cylin- 
drical surface  area  equals  the  spherical  surface  area  (2^r  -.1=4^2  .). 

(S0)cyl     =    iVJal  (^sphere    "     iVsphere 

0  0 
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Figure  20.  Color  photograph  of  liquid  crystals  in  thermal 
field  of  cylindrical  probe  (L/D  =  4).  Current  flux  =  0.414 
amps/cm2 


Figure  21.  Color  photograph  of  liquid  crystals  in  thermal 
field  of  cylindrical  probe  (L/D  =  10).  Current  flux  =  0.414 
amps/cm? 
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Figure  22.  Color  photograph  of  liquid  crystals  in  thermal 
field  of  cylindrical  probe  (L/D  =  20).  Current  flux  =  0.414 

amps/cm? 


Figure  23.  Color  photograph  of  liquid  crystals  in  thermal 
field  of  spherical  probe 
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Figure  24 
Comparison  of  cylindrical  data  (L/D  =  4)  and  cylindrical 
theory  with  spherical  theory  and  scaled  cylindrical  data 
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v  oy  cyl     =    cyl 
/c  \    fr^ 

v  o'sphere      sphere 

but  r2  .    =  r  -,L/2  from  the  equivalence  of  surface  areas.  There- 
sphere    cyl  '  ^ 

fore, 

^'sphere  =  4(Go)cyl 
Similarly, 

r  =^L  T  =  Zal 

sphere     2  sphere      4 

As  can  be  seen,  the  experimental  data  is  more  closely  approximated  by 
the  spherical  theory. 

The  cylindrical  probe  with  a  length  to  diameter  ratio  of  20  most 
closely  approximated  the  infinite  probe  of  the  cylindrical  theory. 
Use  of  the  correction  factor  to  account  for  current  path  changes  near 
the  air-agar  interface,  tends  to  bring  theory  and  experiment  into 
better  agreement. 

B.   RECOMMENDATIONS 

Further  work,  both  experimental  and  theoretical,  is  needed  for  the 
cylindrical  geometry.  Of  particular  interest  are  the  cylindrical  probes 
with  length  to  diameter  ratios  of  five  or  less  since  actual  probes  will 
be  in  this  range.  Experimental  investigations  into  the  effect  of  the 
air-agar  interface  with  varying  probe  depths  would  demonstrate  the 
analogous  effects  of  lesion  creation  close  to  the  skull;  both  the  air- 
agar  interface  and  the  skull  are  essentially  insulators. 

Liquid  crystals  placed  perpendicular  to  the  axis  of  the  cylindrical 
probe  and  running  from  the  air-agar  surface  to  the  bottom  of  the  agar 
could  effectively  map  the  distorted  thermal  region  due  to  the  insulated 
air-agar  surface. 
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Studies  of  various  tissues  should  be  conducted  to  fully  catalogue 
tissue  electrical  properties.  Such  a  catalogue  would  allow  the  surgeon 
to  theoretically  determine  current  strength  and  duration  for  a  particu- 
lar lesion  size  prior  to  the  actual  surgery. 

Some  refinement  in  the  spherical  experimentation  would  reduce 
uncertainty.  The  major  uncertainties  were  in  the  thermophysical  prop- 
erties of  the  agar. 

A  method  of  retaining  the  liquid  crystals  in  the  plane  of  a  diameter 
of  the  probe  is  required.  This  is  particularly  difficult  when  using  the 
cylindrical  probe.  Placing  the  substrate  in  tension  in  its  holder  only 
succeeds  temporarily.  The  substrate  of  mylar  eventually  stretched  and 
sagged.  Without  having  the  liquid  crystals  in  the  diametral  plane,  true 
distances  cannot  be  measured  for  scaling  with  respect  to  the  probe  radius 

A  better  potentiometer  is  required  for  controlling  the  RFG-2AV  gen- 
erator output.  The  control  provided  by  the  manufacturer  had  excessive 
play.  This  was  partially  overcome  by  careful  selection  of  meter  scales 
during  experimentation. 

Use  of  a  35mm  camera  with  a  through  the  lense  light  meter  system  is 
suggested  for  taking  data.  Most  data  for  this  study  was  recorded  using 
the  industrial  Polaroid  film.  When  used  under  artificial  light,  the 
film  is  wery   slow  (ASA  24);  colors  are  not  well  defined  which  is  criti- 
cal when  using  liquid  crystals.  The  net  result  with  the  color  Polaroid 
is  a  poor  reproduction  of  the  liquid  crystal  colors.  Using  35mm  black 
and  white  or  tungsten  corrected  color  film  gives  highly  satisfactory 
results  when  used  with  the  proper  light  meter  settings.  Cameras  using 
Polaroid  film  are  often  moved  during  film  insertion  and  removal.  The 
slightest  movement  results  in  an  unfocused  picture  at  the  short  distances 
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normally  used  in  photographing  the  probe  area.  With  a  35mm  camera, 
there  is  a  minimum  of  movement  during  the  shutter  release  and  film 
winding  processes.  The  actual  time  for  inserting  and  removing  the 
Polaroid  film  limits  the  data  that  can  be  taken;  this  is  most  true  at 
small  times.  Using  a  35mm  camera,  with  the  fast  film  advance  time, 
allows  more  frequent  data  recording. 

A  study  should  be  conducted  to  determine  the  feasibility  of  mix- 
ing more  than  three  liquid  crystals  without  a  loss  of  color  brilliance, 
The  required  separation  of  temperature  ranges  should  be  determined 
since  the  closer  the  temperature  ranges,  the  more  crystals  could 
theoretically  be  mixed.  This  would  yield  more  data  for  each  experi- 
mental trial.  The  temperature  ranges,  however,  must  be  far  enough 
apart  so  that  each  color  band  is  displayed  clearly. 
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APPENDIX  A:  RESISTIVITY 

BACKGROUND 

Initially,  thermal  and  electrical  properties  of  the  99.7  percent 
water/0.3  percent  agar  solution  were  unknown.  It  was  decided  to  inves- 
tigate the  electrical  resistivity  of  the  solution  since  thermal  proper- 
ties could  be  correlated  with  the  solution's  water  content  [10].  A 
test  cell,  pictured  in  Figure  25,  was  designed  and  manufactured  for 
use  in  the  resistivity  investigation.  To  prevent  a  standing  wave  in 
the  resistivity  test  cell,  the  test  cell  length  could  not  be  an  integral 
or  a  quarter  wavelength  of  the  500  kilohertz  signal  used  in  the  experi- 
ments. Using  the  relationship  between  frequency  (f),  wavelength  (a), 
and  the  speed  of  light  in  air  related  to  the  speed  of  light  in  water 
(n), 

3  x  108 
x  =   J  x  'u    meters 

the  wavelength  was  determined  and  the  test  cell  designed  to  prevent 
standing  waves. 

A  plexiglass  tube,-  13  cm.  in  length  with  2.5  cm.  ID,  was  used  as 
the  test  cell  due  to  the  insulating  properties  of  plexiglass  at  the 
power  levels  used.  Copper  plugs  were  machined  to  fit  tightly  into  the 
cylinder  with  the  area  of  the  plug  ends  known.  A  BLH  two- inch, 
constantan-chromel  needle  thermocouple  was  introduced  through  the  wall 
of  the  plexiglass  cylinder  midway  between  the  copper  ends.  To  prevent 
leakage  of  air  into  the  container,  a  seal  of  neoprene  rubber  was  provided 
between  the  test  cell  and  the  thermocouple  shaft.  During  the  experimental 
trials,  the  test  cell  was  filled  with  hot  agar  solution. 
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By  carefully  capping  the  test  cell,  no  air  was  trapped  inside. 
The  capped  cylinder  was  placed  in  the  generator  circuit  shown  in 
Figure  26. 

Voltage  was  applied  across  the  circuit  for  one  minute  when  col- 
lecting data.  One  minute  was  chosen  since  no  change  in  voltage  drop 
occurred  after  this  time.  Determination  of  resistivity  was  accom- 
plished by  recording  voltages  across  the  circuit  and  the  ten  ohm 
resistor  along  with  the  distance  between  the  surfaces  of  the  copper 
plugs.  The  voltage  drop  across  the  resistor  yielded  circuit  current 
from  Kirchoff's  Law: 

I  = 


R10 

The  resistance  of  the  cylinder  was  determined  by  dividing  this  current 
into  the  voltage  drop  across  the  cylinder.  Resistivity  is  related  to 
resistance  by  the  area  to  length  ratio: 

-   A    R 


M      L 
It  was  thus  possible  to  calculate  resistivity. 

At  the  same  time  the  voltage  readings  were  recorded,  the  needle 
thermocouple  was  read  via  the  millivolt  potentiometer.  Data  were  taken 
as  the  test  cell  cooled  and  again  as  the  cell  was  heated  in  a  controll- 
able oven  (Figure  27).  Early  experimental  trials  showed  a  wide  variation 
in  agar  resistivity  (Figure  28).  To  reduce  all  plots  to  a  common  scale, 
two  nondimensional  quantities  were  introduced: 

T  -  "L 

Temperature 


T 
o 

p  =  -e—  Resistivity 

po 


60 


MILLIVOLT  POTENTIOMETER 
ICE  BATH 


RESISTIVITY  TEST  CELL 


RADIO  FREQUENCY 
GENERATOR 


Voltmeter  V-.  measured  voltage  drop  across  the  ten  ohm  resistor 
from  which  circuit  current  v/as  determined. 

Voltmeter  V?  measured  circuit  voltage  which  was  used  in  conjunction 
with  circuit  current  to  determine  resistivity  test  cell  resistance. 


Resistivity  Test  Cell  Circuitry 
Figure  26 
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Figure  27.  The  resistivity  test  cell  in  the  oven 
for  the  heating  phase  of  resistivity  data  collection 
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Figure  28 
Variation  in  resistivity  values 
as  a  function  of  temperature 
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where  the  subscript  indicates  the  initial  conditions  of  the  experimental 
trial.  In  this  way,  the  plots  had  common  points  of  origin. 

A  linear  curve  was  fit  to  each  set  of  data  points.  In  dimensional 
form, 

P=Po  +  Y  (T.-To) 

where  y  was  the  slope  of  the  curve.  Nondimensionalizing  yielded 

p  =  1 •+  me 

yT 
where  m  =  — -  ,  the  nondimensional  slope.  The  values  of  m  were 

po 

generated  by  computer  program  III  (Appendix  E).  Figure  29  is  a  sample 
of  the  nondimensional  curves  plotted  by  this  program. 

As  a  result  of  this  portion  of  the  investigation,  the  great  change 
in  resistivity  over  the  temperature  range  of  interest  (15-60°C)  was 
demonstrated.  It  was  determined  that  the  programs  written  by  Gengler  [5] 
required  modification  to  reflect  the  changing  resistivity.  This  was 
accomplished  by  incorporating  the  slope  of  the  linear  curve  from  the 
resistivity  data  into  computer  programs  I  and  II  (Appendix  E). 
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APPENDIX  B:  LIQUID  CRYSTALS 

BACKGROUND 

Matter  can  exist  in  solid,  liquid,  or  gaseous  state  depending  upon 
the  temperature  and  pressure  of  the  surroundings.  In  the  solid  state, 
most  compounds  are  composed  of  an  ordered  array  of  molecules  in  a 
crystal  structure.  When  heated,  the  crystal  begins  to  deform  as  the 
molecules  absorb  energy  and  leave  the  "ground  state."  The  temperature 
at  which  the  deformation  begins  is  known  as  the  melting  point  [17]; 
it  is  normally  a  well-defined  characteristic  of  the  compound. 

In  1888  Reinitzer  noticed  that  not  all  compounds  changed  directly 
from  solid  state  to  an  isotropic  liquid  [18].  The  intermediate  phase, 
a  "mesophase,"  exhibited  characteristics  of  both  crystalline  structure 
and  isotropic  liquid.  Thus,  the  name  of  these  substances  was  coined 
as  liquid  crystals  [18].  Lehmann  continued  Reinitzer's  work  in  liquid 
crystals.  He  classified  the  crystals  into  three  categories  based  on 
molecular  orientation;  the  categories  are  smectic,  nematic,  and 
cholesteric  [17].  Only  the  cholesteric  crystals  will  be  discussed  as 
it  is  the  only  type  utilized  in  this  report.  Interested  readers  may 
find  information  on  smectic  and  nematic  crystals  in  the  literature 
[16,  17,  18,  19,  20]. 

The  cholesteric  mesophase  has  the  molecules' long  axis  parallel  and 
coplanar.  A  model  for  this  structure  is  a  deck  of  cards,  each  card 
with  one  corner  bent  up  such  that  a  slight  rotation  of  the  entire 
stack  results  [18].  The  average  displacement  of  adjacent  layers  is 
fifteen  minutes  of  arc  per  layer  [21]  with  each  layer  on  the  order  of 
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1000  A  thick  [22].  The  rotation  is  due  to  the  stretched  out  non- 
symmetric  shape  of  the  crystal  molecule  [23]. 

The  helical  rotation  of  the  cholesteric  liquid  crystal  gives  rise 
to  its  optical  properties.  When  transmitting  polarized  light,  the 
helical  layers  cause  the  light  to  be  rotated  to  the  left  with  rotations 
of  fifty  degrees  per  millimeter.  This  property  is  known  as  being 
optically  active.  By  way  of  comparison,  quartz,  one  of  the  most 
optically  active  materials  in  nature,  has  a  rotation  of  twenty  degrees 
per  millimeter  for  polarized  light.  Another  property  of  the  cholesteric 
liquid  crystals  is  circular  dichroism  which  is  a  function  of  temperature, 
material,  and  the  angle  of  the  incident  light.  The  incident  beam  of 
unpolarized  white  light  is  split  into  two  components.  One  component 
is  rotated  clockwise  while  the  other  component  is  rotated  counter- 
clockwise. The  result  is  that  one  component  is  reflected  while  the 
other  is  transmitted.  To  absorb  the  transmitted  beam,  a  black  background 
is  required  behind  the  liquid  crystals.  Circular  dichroism  is  the  char- 
acteristic that  gives  cholesteric  liquid  crystals  their  iridescent 
colors.  This  characteristic  should  not  be  confused  with  the  birefringence 
or  double  defraction  of  white  light  by  the  cholesteric  liquid  crystals 
(Figure  30).  Birefringence  is  often  used  as  a  test  for  the  cholesteric 
mesophase  [21].  In  the  material  the  incident  beam  is  divided  into  two 
polarized  components  which  are  vibrating  at  right  angles  to  each  other 
and  are  travelling  at  different  speeds.  They  are,  therefore,  refracted 
at  different  angles  and  emerge  as  two  beams  of  perpendicularly  polarized 
light.  The  rotary  power  of  the  birefringent  cholesteric  mesophase  is 
as  much  as  1800  degrees  per  centimeter  when  light  is  transmitted  along 
the  optical  axis  normal  to  a  thin  sheet  [23,  24].  The  thickness  of  the 
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Birefringence:  In  the  material,  the  incident  beam  of  light  is  divided 
into  two  polarized  components  vibrating  at  right  angles  to  each  other 
and  traveling  at  different  speeds.  They  are,  therefore,  refracted  at 
different  angles  and  emerge  as  two  beams  of  perpendicularly  polarized 
light 


Circular  dichroism:  The  incident  beam  of  unpolarized  light  is  split 
into  two  components,  one  rotating  clockwise  and  the  other  rotating 
counterclockwise.  The  result  is  that  one  is  reflected  while  the  other 
is  transmitted. 


Figure  30 
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cholesteric  film  does  not  alter  the  predominant  wavelength  of  the 

reflected  light  as  it  is  rotated;  however,  the  angle  of  incidence 

and  the  angle  of  observation  both  effect  the  wavelength  of  the  reflected 

light. 

The  most  striking  changes  in  the  cholesteric  mesophase  are  those  due 
to  temperature  changes.  This  mesophase  is  made  temperature  sensitive 
by  the  narrow  band  of  optical  reflection  which  is  associated  with  a 
change  in  crystal  structure.  Cholesteric  liquid  crystals  are  also 
sensitive  to  organic  compounds,  shear,  and  electromagnetic  radiation 
at  high  frequencies  and  power  levels.  Depending  on  composition,  chol- 
esteric liquid  crystals  span  a  range  from  -20  to  250  degrees  C  [25]. 
Any  specific  compound  may  have  its  range  narrowed  to  as  little  as  one 
degree  or  increased  to  six  degrees  by  varying  the  compound  ingredients 
[26].  If  heated  far  beyond  its  temperature  range  and  allowed  to  cool, 
the  transition  colors  occur  at  lower  temperatures,  similar  to  a  magnet 
being  heated  above  the  Curie  temperature  [21].   If  not  thermally  shocked 
as  just  described,  the  transition  colors  remain  invariant  with  tempera- 
ture; the  material  may  be  repeatedly  recycled.  Not  all  cholesteric 
liquid  crystals  react  with  the  same  color  changes;  some  display  all  colors 
of  the  spectrum  starting  with  red  to  yellow  to  green  to  blue,  while  others 
may  change  directly  from  red  to  blue. 

Cholesteric  liquid  crystals  will  deteriorate  in  several  days  if  not 
insulated  from  dust  and  organic  materials.  National  Cash  Register 
Company  has  been  successful  with  encapsulating  the  cholesteric  liquid 
crystals  in  spheres  5  to  30  microns  in  diameter  [23].  The  result  of 
the  microencapsulation  process  is  a  long  life  with  little  deterioration 
from  contamination.  The  NCR  process  reduces  the  angle  dependence  of  the 
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crystal  as  far  as  color  is  concerned  although  resolution  is  degraded 
and  response  time  is  increased  approximately  ten  percent  [22].  The 
microencapsulated  crystals  are  used  in  industry  since  they  are  not 
readily  contaminated;  they  were  used  in  this  experiment. 

Liquid  crystals  have  many  engineering  applications.  In  temperature 
measurement,  cholesteric  liquid  crystals  are  primarily  used.  Guidelines 
for  use  of  liquid  crystals  in  thermal  mapping  may  be  summarized  as 
follows: 

a)  The  heat  capacity  of  the  object  should  be  greater  than  that 
of  the  liquid  crystals; 

b)  The  size  of  the  object  must  be  great  enough  to  resolve  the 
colors  of  the  liquid  crystals; 

c)  The  rate  of  temperature  change  must  be  slow  enough  for  the 
liquid  crystal  to  follow.   (A  maximum  time  of  about  0.2  seconds.); 

d)  The  temperature  range  under  study  must  be  within  the  range  of  the 
liquid  crystal ; 

e)  The  surface  of  the  object  must  be  capable  of  holding  liquid 
crystals  unless  an  adhesive  backing  is  used; 

f)  The  surface  must  be  black  to  display  the  liquid  crystals; 

g)  The  temperature  pattern  must  be  correlated  with  another 
material  property  such  as  voids  present  in  nondestructive 
testing  or  power  dissipated  per  unit  area. 
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APPLICATION  AND  CALIBRATION 

Reference  27  lists  many  techniques  for  application  of  liquid 
crystals;  in  this  investigation,  a  commercial  artist's  air  brush  was 
utilized.  Three  crystals,  designated  R-27,  R-37,  and  R-53,  were  mixed 
equally.  After  standing  overnight,  the  cholesteric  substances  tended 
to  separate  from  excess  water.  The  water  was  carefully  removed  leaving 
a  more  concentrated  liquid  crystal  which  resulted  in  a  more  brilliant 
color  display.  The  concentrated  mixture  was  sprayed  onto  the  transparent 
substrate  in  three  thin  coats,  allowing  adequate  drying  time  for  each 
coat.  To  prevent  any  absorption  of  moisture  from  the  environment,  dry- 
ing was  conducted  under  a  40  watt  lamp.  After  the  three  coats  of  liquid 
crystals  were  dried,  two  coats  of  black  paint  were  applied.  Testors 
Flat  Black  Spray  Enamel  was  used  since  it  contained  no  chemical  which 
might  have  reacted  with  the  crystals  to  cause  a  shift  in  the  color- 
temperature  relationship.  The  black  paint  was  required  to  absorb  the 
light  transmitted  through  the  liquid  crystal,  as  only  a  small  portion  of 
the  light  was  reflected.  To  provide  a  watertight  system,  a  layer  of 
polyurathane  varnish  was  applied  over  the  black  paint  and  along  each  edge 
of  the  substrate.  Thus,  the  crystals  were  sealed  from  all  contaminants. 

A  second  type  of  liquid  crystal,  designated  ROCHROME,  was  also  used 
during  the  experiments.  Rochrome,  manufactured  by  Hoffman-LaRoche, 
consists  of  liquid  crystals  sandwiched  between  a  layer  of  mylar  and  a 
black  adhesive  backing.  Rochrome  is  cut  to  the  desired  shape  and  applied 
directly  to  the  object  under  study.  Again,  the  edges  were  sealed  to  pre- 
vent water  contamination  and  separation  of  the  Rochrome  from  the  surface. 
The  surfaces  to  which  the  Rochrome  and  the  liquid  crystals  were  applied 
in  this  investigation  were  mylar  and  glass.  The  mylar  was  used  with  the 
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probe  apparatus  while  the  glass  was  used  during  the  calibration 
procedure. 

Calibration  was  conducted  using  the  Rosemount  Engineering  Company 
Variable  Temperature  Bath  (Model  913A)  in  conjunction  with  its 
commutating  bridge  and  platinum  resistance  standard  thermometer. 

National  Cash  Register  liquid  crystals  (R-27,  R-37,  and  R-53) 
were  prepared  by  spraying  them  onto  a  glass  slide,  coating  with 
black  paint,  and  sealing  with  polyurathane  varnish.  Rochrome  was 
prepared  by  sticking  a  small  piece  to  the  slide  and  sealing  the  edges 
with  polyurathane.  The  slides  were  placed  in  the  bath  at  approximately 
one-third  of  the  bath's  depth  and  close  to  the  center  of  the  bath 
diameter.  Effects  of  the  bath  container  edges  were  thus  minimized. 
The  bath  temperature  was  raised  slowly  as  the  slides  were  observed. 
To  minimize  variations  between  observers,  the  same  observer  read  the 
slide  colors  while  an  assistant  operated  the  equipment  during  all 
calibration  procedures. 

As  the  temperature  was  raised,  the  liquid  crystals  reacted,  each 
at  its  own  temperature  level.  When  the  slide  had  turned  one  color, 
i.e.  red  to  green,  the  commutating  bridge  reading  was  recorded;  then 
the  temperature  of  the  bath  was  increased  to  determine  the  next  color 
change.  A  conversion  table  was  used  to  convert  the  bridge  readings  to 
temperature.  The  uncertainty  in  this  temperature-color  calibration 
procedure  was  calculated  to  be  ±0.5°C  which  agrees  with  Reference  27. 
Table  B-l  lists  the  crystals  used  and  their  calibration  temperature- 
color  relationships. 


72 


Table 

B-l 

Crystal 

Red 

Green 

Blue 

R-27 

26.7°C 

27.9°C 

28.9°C 

Rochrome 

29.9°C 

31.5°C 

34.2°C 

R-37 

35.7°C 

37.0°C 

38.5°C 

R-53 

50.5°C 

52.1°C 

53.3°C 
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APPENDIX  C:  FLUX  PLOTS 

Gengler's  theoretical  analysis  of  the  cylindrical  probe  modeled  the 
probe  as  an  infinite  cylinder  with  negligible  end  effects;  the  current 
flow  should  have  been  a  one-dimensional  problem.  With  the  experimental 
cylindrical  probe,  end  effects  were  pronounced,  prompting  an  investiga- 
tion of  current  paths  in  the  probe  vicinity.  The  area  of  interest  was 
within  ten  probe  diameters.  The  problem  was  approached  in  two  phases. 
First,  the  horizontal  plane  was  examined  using  TRUMP  [28,  29],  a  computer 
code,  while  in  the  second  phase,  teledeltus  paper  was  used  to  depict  the 
isopotentials  in  the  vertical  plane. 

TRUMP  solved  the  general  partial  differential  equation  v2e  =  0. 
It  was  postulated  that  the  radio  frequency  potential  steady  state 
condition  was  established  instantaneously  with  power  application. 
A  nodal  network  was  constructed' for  use  with  TRUMP  (Figure  31).  Due 
to  symmetry  of  the  probe-cell  system,  only  one  quadrant  of  the  system 
was  modeled.  Three  cylindrical  probes,  with  length  to  diameter  ratios 
of  4,  10,  and  20  were  used  with  two  sets  of  boundary  conditions.  The 
first  boundary  conditions  approximated  a  small  finite  cylinder  inside 
an  infinite  cylinder  by  the  use  of  two  conductive  sides  of  the  nodal 
boundaries  for  the  outer  cylinder  wall  and  two  insulated  sides  of  the 
nodal  boundaries  for  the  endcaps  of  the  outer  cylinder.  Figures  32, 
33,  and  34  are  plots  of  the  isopotentials  and  lines  of  constant  current 
for  each  probe  length  to  diameter  ratio.  Inspection  of  these  figures 
in  the  area  within  ten  probe  diameters,  shows  the  current  lines  orthogonal 
to  the  probe  surface.  It  was,  therefore,  expected  that  the  current  flow 
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TRUMP  nodal  network 
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and  the  resultant  thermal  problem  would  be  one  dimensional.  For 
comparison,  a  second  set  of  boundary  conditions  surrounded  the  probe 
with  conductive  surfaces.  Figures  35,  36,  and  37  display  the  iso- 
potentials  and  constant  current  lines  for  this  situation. 

The  second  portion  of  current  path  investigation  was  undertaken 
when  experimental  data  did  not  agree  with  the  theory.  As  stated  earlier, 
in  this  second  phase,  isopotentials  in  the  vertical  plane  were  examined 
using  the  same  sets  of  boundary  conditions  as  in  the  first  phase. 

The  probe-cell  system  was  scaled  on  teledeltus  paper.  Boundary 
conditions  were  simulated  using  conductive  silver  paint  and  metallic 
tape.  Figures  38  and  39  are  plots  of  the  isopotentials  and  constant 
current  lines  under  the  two  sets  of  boundary  conditions.   It  is  seen 
that  the  insulated  agar  surface  caused  current  to  bend  towards  the 
conductive  boundary.  Since  heating  is  a  function  of  current  per  unit 
area,  the  higher  current  density  resulted  in  higher  heating  effects 
in  the  experimental  model.  The  "higher  heating  effects  were  the  major 
cause  for  the  difference  between  theory  and  experiment  with  the 
cylindrical  probe. 
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<j>'  =  0  (Free  surface  of  agar-water  mixture) 


Figure  38 

Teledeltus  plot  of 
isopotentials  and  constant 
current  lines  in  a  plane 
perpendicular  to  the  cyl- 
indrical probe  axis. 

Free  surface  of  agar- 
water  mixture  acting  as  an 
insulator. 


(Free  surface  of  agar-water  mixture) 

(J)  =  0 


Figure  39 
Teledeltus  plot  of 
isopotentials  and  con- 
stant current  lines  in 
a  plane  perpendicular 
to  the  cylindrical 
probe  axis. 

«j>'-G 
Free  surface  of  agar- 
water  mixture  acting 
as  conductive  boundary. 


4>  =  0 
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APPENDIX  D:  SAMPLE  UNCERTAINTY  ANALYSIS 

Reference  30  categorizes  error  in  experimental  investigation  as 
fixed  error,  variable  error,  and  mistakes.  Fixed  error  is  defined 
as  a  constant  error  such  as  equipment  consistently  reading  low. 
Variable  error  is  caused  by  variations  in  sensors,  an  example  of 
which  is  having  two  different  observers  taking  readings  in  separate 
experimental  trials.  Finally,  mistakes  may  occur  during  experiments 
by  reading  the  wrong  scale,  etc.  If  it  is  assumed  that  the  error, 
i.e.  uncertainty,  of  each  primary  variable  has  the  same  probability 
of  occurrence,  it  is  then  possible  to  calculate  the  overall  uncer- 
tainty in  each  value.  As  an  example,  define  P,  the  result  of  some 
experiment  as , 

If  to  is  the  uncertainty  in  P,  and  w-. ,  w«,  u^,  ...w  are  the 
uncertainty  in  each  primary  independent  variable,  then  the  relative 
uncertainty  in  P  is  defined  as 


J*e. 


In  the  experimental  trials  on  the  spherical  probe,  the  quantities 
calculated  were 


G 

o 


i2po  v2po 

16^r  2k  T         =         16^r%  T  O 
oo  o  o    10 
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T  -  T 


T   = 


at 


R  = 


Deriving  the  expressions  for  the  relative  uncertainty  in  each 
quantity 


U) 


G      , 
o 


w. 


'2u).  ^2 


T  -  T 


p  * 

l-£l  + 


(0T      2 
o 


"::r  •  r.  • 


fuT\2 


f2(0D  v2 


w 


1/2 


1/2 


V2 


(0 


V 


1/2 


For  all  independent  variables  except  resistivity  p   it  was  possible 
to  estimate  the  uncertainty  either  directly  from  observing  equipment 
(voltage,  current)  or  from  available  tables  (a,  k).  For  resistivity 
it  was  necessary  to  calculate  uncertainty  using 


■•  ■  m  ■  ^  '" 


fr)ff) 


Following  the  example  of  Equation  D-l , 
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u 


fO)i 


cct 
Vcct  -  VR 


'10 


Vcct  -VRI 


cct 

r 

cct 


10 


1/2 


Table  D-l  lists  values  of  uncertainty  for  each  independent  variable 
and  the  value  of  the  variable  used  in  the  calculations.  The  following 
values  were  calculated  for  uncertainty: 


CO 


0.036   (3.6%) 


U)r 


0.1217   (12.2%) 


U), 


0.019   (1.9%) 


R 


0.029   (2.9%) 


CO 


0.102   (10.2%) 


To  obtain  uncertainty  in  the  quantity  e/G  the  following  relationship 


was  used 


0J 


e/G 


e/G 


el   ./  o1 


Using  the  values  above, 


1/2 


CO 


e/G 


e/G 


0.123   (12.3%) 
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As  may  be  seen  from  the  last  column  of  Table  D-l,  the  largest 
uncertainties  arise  from  the  thermophysical  properties,  i.e.  thermal 
conductivity  (k),  thermal  diffusivity  (a),  and  resistivity  (po). 


Table  D-l 

Quantity   i 

Uncertainty  wi 

Percent  Uncertainty 

V 
vcct 

=  80.0  volts 

0.0025  volts 

0 

VR 

=  0.118  volts 

0.0025  volts 

0.02 

D 

=  2.540  cm 

0.025  cm 

0 

L 

=  11.95  cm 

0.025  cm 

0 

R10 

=10.0  ohms 

0.3  ohms 

0.03 

V 

=  0.60  volts 

0.0025  volts 

0 

r 

0 

=  0.241  cm 

0.001  cm 

0 

T 

0 

=  23.0  °C 

0.25  °C 

0.01 

T 

=  37.0  °C 

0.10  °C 

0 

r 
°P 

rp 

=  0.380  cm 
=  0.90  cm 

0.01  cm 
0.01  cm 

0.02 
0.01 

t 

=  30.0  sec 

0.025  sec 

0 

k 

=  0.628  W/m°K 

0.0628  W/m°K 

10.0 

a 

=  0.001512  m2/sec 

0.0001512  m2/sec 

10.0 

po 

=  3250.0  ohm-cm 

117.0  ohm-cm 

3.60 
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APPENDIX  E 

C  COMPUTER    PROGRAM     I 

C  PROGRAM    FOR    THE     INSULATED    SPHERE. 

C 

c 

C      MODIFIED  TO  REFLECT  CHANGING  RESISTIVITY 
C 

c 
c 

C  TET    REPRESENTS    THE    NONDIMENS ION AL    TEMPERATURE. 

C  BETA    IS    THE    BLOOD    FLOW    TERM. 

C  G     IS    THE    NONDIMENSIONAL    POWER    TERM. 

C  DR    IS    DELTA    RADIUS. 

C  EM    REPRESENTS    THE    M    IN     THE     STABILITY    CRITERIA. 

C  DELTAU    REPRESENTS    NOND IN  EMS  ION A L    TIME    CHANGE. 

DIMENSION    TFT (50 0.2) ,TT( 500,2) , RH0*500) 

ITIMF=0 

R E AD ( 5. 100) BET A, DR. GO. AM, LI, I H 
100    FORMAT* 3F10.4,F10.6,2I10) 

WRITE(6. 1 50 ) BET  A, DR, GO, AM,  LI , IH 
150    FORMAT( • 1« .6X,« 3ETA  =  • , F4 . 2, / , 7X, ' DR  =  • , F6. 4, /, 7X, • G 
AO' =  • ,F6. 1,/,7X, • SMALL  M  = • , Fl 0 . 6 , / , 7X , • L I  =  «,I3,/t7X 
B.'HI  =  '.1.4.//) 

TET*1  .  1  )=0.0 

WRITE(  6,2000) 
2000       FORMAT (• l',T4. 'TEMPERATURE    DISTRIBUTION    AT    SPECIFIED    T 
AIME    AND    RADIAL    LOC AT  IONS ■  , / / , 9 X , • T  IME     ( T AU  )  ■ ,  IX ,  • LOCAT 
BION    (R)  ',  3'X.  'TEMPERATURE    (  THETA  ) /PCWER     (GO)') 

S=AM*DR**2#GQ 

E  M=  1.  0  /  <  2-  0+ BETA^  DR**2-S ) 

DELTAU- EM^  DR*  ■  Z 

DO    10     K  =  1..IH 
10  TFT*K,l)=Q.O 

50  DO    200     I=LI. IH 

TET(IH+1.1)=TET{ IH-1,1 ) 

TFT*LI-1, L)=TET*  LI  +  1 ,1) 

X=FLGAT ( I) 
200  TET  (I  ,2)=EM-'  (  ( 1 .  -1  .  /X  )  *TET  (  1-1,1  )  +(  1  .  +1  .  /  X  )#T  ET  (  1  +  1,  1) 

R)+TET(  I  «  1  )■•.•(  l.-2.*EM+(  EM*DR**2* AM*GO> / <  X*DR)**4)+  (  GC:- E 
CM*DR**2  )/4X*.  DR  )**4     ' 

I TIMF-TTIMF+1 

TIME= FLOAT* I  TIME) 

TAU-TIMEH-DELTAU 

IFUTIME.GT.1000)     GO    TO    90 

IF*  IT  IMF.GT.100)     GO    TO    70 

•ITEST  =  ITIME/20*20 

GO    TO     8  0 
70  ITEST= IT IME/ 50*50 

GO    TO    80 
90  ITFST=ITIME/100*100 

80  J=0 

DO    60    J  =  L1  .  IH 

Y=FLOAT* J) 

R=Y*OR 

IF* * ITIME-ITEST) .NE.O)     GO    TG    60 

TT* J.2)=TET*  J, 2)  /GO 

WRITE(6.220)     T A J , R, TT *  J  ,  2  ) 
220  FORMAT*'     «..  Fl  5.  4  .  7X,  F6  .3  ,  5X  ,  El  5  .7.) 

60  TET*J. 1)=TET( J, 2) 

IFUTIME.GT.4000  )STOP 

GO    TO     5  0 
END 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


100 
150 

2000 


10 
50 


200 


COMPUTER  PROGRAM  II 

PROGRAM  FOR  THE  INSULATED  CYLINDER 

MODIFIED  FOR  VARIABLE  RESISTIVITY 


70 

90 

80 


220 
60 


TFT  RE 
BETA  I 
G  IS  T 
DR  IS 
EM  REP 
DELTAU 
D  I  MENS 
ITIMF= 
READ(5 
FORMAT 
WRITc( 
FORMAT 
AO  =  •  . 
BHI  =  ■ 
WR I TE ( 
FORMAT 
A  I  ME  AN 
BION  ( R 
FM=1.0 
DFLTAU 
DO  10 
TFT(K. 
DO  200 
TET(  IH 
TET(LI 
X=I 

TETU, 
ATETt  1  + 
BET(  I,  1 
IT  IMF* 
TAU=IT 
I F(  IT  1 
I F(  IT  T 
ITEST= 
GO    TO 
ITEST= 
GO    TO 
ITEST= 
J=0 
DO    60 
R=J4DR 
IF( (IT 
TT( J, 2 
WRITE( 
FORMAT 
TFT( J, 
IF(  IT  I 
GO    TO 
END 


PRES 
S  TH 
HE  N 
DELT 
RE  St 

RFP 
ION 
0 

,  100 
(  3F1 
6,15 
(«  1«. 
F6.3 

14* 
2  0 

1  ie: 

RA 

•  ,3 

(2. 

=  EM.J 

K  =  l, 

1  )  =  0 

I  =  L 

+  1.1 

-l.L 


)BET 

0.4, 

OJBE 

,6X. 

,/.7 

//) 

00) 

,T4. 

DIAL 

X.  "T 

0-3E 

DR*-* 

IH 

-0 

I.IH 

)  =  TF 

)  =  TE 


NONDIMENSIONAL    TEMPERATURE. 
FLOW    TERM. 
IONAL    PDWER    TERM. 

M    IN    THE     STABILITY    CRITERIA. 
NONDIMENS IONAL    TIME    CHANGE. 
2)  .TT<500,2> 

A, DR. GO. AM, LI. IH 

F10.6.2I  10) 

TA.DR.GO, AM.LI , IH 

•BETA    =     • ,F4.2,/,7X,»DR    =     « , F6. 4. / , 7X, • G 

X, 'SMALL    M    =     • .F10.6,/.7X, • LI    =    ',I3»/,' 


ENTS    THE 
E    BLOOD 
ONDIMENS 
A    RADIUS 
NTS    THE 
RESENTS 
TFT( 500. 


•TEMPERATURE    DISTRIBUTION    AT    SPECIFIED    T 
LOCATIONS' ,//,9X. "TIME    (TAU)  ', IX, 'LOCAT 
EMPERATURE     ( THET A ) /POWER     (GO)') 
TA*DR**2-AM*G0) 
2 


2)  =  EM::  (  (  l.-l 

1 ,1 ) >+(  l.-BE 

)+(GO--EMJ  /X* 

ITIME+1 

I ME&DELTAU 

ME. ST. 1000) 

ME.GT.100J    G 

I  TIME/20*20 

80 

ITIME/50*50 

SO 

ITIME/IOCX-'IO 

J=LI.IH 

IME-ITEST  )  .N 
)=TET< J. 2) /G 
6,220)TAU.R. 
( •  »  ,F15.4,7 
1)=TET( J,2) 
ME.GT.4000JS 
50 


T(IH-l.l) 
T(LI+1, 1 ) 

,/<2.*X) )*TET( 1-1,1) +(1.+1./ <2.*X))* 

TA*  EM*DR#*2-2.<  EM+( AM*GG*  EMI/X**2)*T 
*2 


GO    TO    90 
0    TO    70 


E.O  )     GO    TO    60 

Q 

TTU.  2) 
X.F6.3.5X,E15.7i 

TOP 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


1 

2 

3 

11 

12 

7 
4 

6 
15 

8 
99 

13 
101 


104 
100 


PROGRAM 
PROGRAM 
TEMPERA 
THE    CUR 
THESE    C 
LSQPL2 
FOR       AL 
M=NUMPT 
TEMP=MO 
RES=RES 
OX    REPR 
DF2     IS 
IMPLICI 
REAL  1*8 
REAL?8 
REALS8 
REAL    LA 
REAL    LA 
D I  MEN  SI 
ADB(IOO) 
B( 100)  .1 
READ( 5. 
FORMAT ( 

no   ioo 

READ( 5, 
FORMAT ( 
WR I TE ( 6 
FORMAT ( 
AHE  MAXI 
BLE  IS  • 
CIS  SCAL 
RFAD( 5, 
FORMAT( 
READ(5. 
FORMAT( 
W  R  I T  t  (  6 
FORMAT ( 
DTG=TO 
DO  7  1  = 
ORES! I) 
DTHET( I 
DW  I  (  I  I  = 
WRITE  (6 
FORMAT ( 
BLUES  .  / 
C  RHO-P 
WRITFC  6 
FORMAT ( 
CALL  LS 
00  15  I 
B(I)=DS 
DO  8  1  = 
Y  (  I  )=DY 
THETM I 
RES(I )  = 
IFIJ.EO 
CALL  OR 
A,  0.  8,  11 
CALL  DR 
AO.8,11. 
WRITE(6 
FORMAT ( 
GQ  TO  1 
CALL  DR 
A. 0,8. 11 
CALL  DR 
A0.8.11. 
WRITE (to 
FORMAT ( 
CO NT  I NU 
STOP 
END 


NUMB 
TO  F 
TURE 
VFS  T 
URVES 
MUST 
L  DRA 
S,  TH 
NDIME 
ISTIV 
ESENT 
RES'IS 
T  REA 
TITLE 
ITITL 
LABEL 
6EL2/ 
3FL1/ 
ON  DX 
tOSBC 
TITLE 
9)  (  IT 
A8) 
J  =  l  ,2 

I  )  M  .  K 
2110. 
,  2  )  M  , 
30  HI 
MUM  0 
F10.1 
E  IS 
3}  (DX 
2F10. 

II  )T0 
F5»  II 
.  12)7 


ER  I 
IT  A 
DATA 
0  DA 
AND 
3EGI 
W  TH 
E  NU 
NSIO 
ITY 
S  NO 
TIVI 
L*8( 
(  10) 
E 
/• 

4HHE 
4HCG 
(  100 
100) 
(12) 
ITLE 


M.  XC 

2F10 

KM.Y 

THE 

EGRE 

,2H 

.F10 

(I)  , 

1) 


I  I 

ND 

f  A 
T 
N 

MB 
NA 
(Y 
NO 
TY 
D) 
/l 


PLOT    CURVES     FOR    RESISTIVITY    AND 
SUBROUTINE    LS^PL2    IS    CALLED    TO    FIT 

.       SUBROUTINE     DRAW    IS    USED    TO    PLOT 

HE    DATA    POINTS.       ALL     ARGUMENTS     FOR 

WITH    D    EXCEPT    M    AND    KM. 

FOLLOWING    NOTATION    IS     USED 

ER    OF    POINTS 

L  TEMPERATURE  (X) 

) 

IMENSIONAL  TEMPERATURE THETA. 


0* '     ■  / 


•/ 


AT/ 

GL/ 

)  .DF2( 100)  ,DWI (100) ,DY(100)  ,DELY(100  ), 

, THET( 100), TEMP ( 100) iRESi 100), B( 100) ,Y 

,DTHET(100) , ORES (100). THETA (100) 

( I )  ,1  =  1,12) 


AL,YCAL 
.1) 

CAL.XCAL 

NUMBER  OF  DATA  POINTS  IS  ,I10,/,30H0  T 

E    OF     FIT     IS     ,I10,/,21H0    THE    Y-AXIS    SCA 

0    IN    DEGREES    PER    INCH     ,/,21H0    THE    X-AX 

.1,20H0  IM  GH^l-CM  PER  INCH  ) 

DF2( I) , 1=1 tM) 


u 
THE 


INITIAL  TEMPERATURE  (TO)  W AS, F5. 1 , 9H0I 


l.M 
=  DF2( 
)  =  (DX 
1.000 
.4) 

14X.4 
,9X,5 
R  1  M  E 
.  6  )  (  0 
12X.F 
0PL2( 
=  1.KM 
(I)     • 
l.M 
(I  ) 
)=DTH 
0RES( 
.2)     G 
AWtiM. 
.  1..LA 
AWiM. 
l.LAS 
.  13) 
8  HO    L 
CO 

AW(M. 
•  1-LA 
AW(M. 
l.LAS 
.1  04) 
8H0  L 
E 


I  I/O 

(  I  )- 

0000 

SHOO 

1H    T 
./,9 

X(  I  ) 
6  .i  , 

M,KM 


ET(I 
I  ) 

G    TO 

THET 

ST) 

THET 

T) 

LAST 

AST 


F211) 
DTO)"/DTO 


IMENSIONAL     VALUES  NCNDI MENSI ONAL    VA 

EMPERATURE       RESISTIVITY  THETA 

X,24H     (DEGREES     C)  (OHM-CM)     ,//) 

,DF2(  I  ) ,DTHET( I  )  .ORES { I  )  , I =1 , M) 
7X,F6.1,8X, F8.o, 5X.F8.6) 
,D THET, ORE S.DWI , DY , DEL Y , DB, CSB ,T ITLE ) 


) 

101 
A, RES, 1,5, LAB  EL. IT ITLE, XCAL , YCAL, 0, 0, 0 

A. Y, 2.0. LAB  EL 2, IT ITLE, XCAL, YCAL, 0,0,0, 


=    ,13) 

A, RES, 2, 4,  LABEL, I  TITLE, XCAL,  YCAL, 0,0,0 

A, Y, 3,0, LAB  EL  1,1  TITLE, XCAL, YCAL, 0,0,0, 


THET 

ST) 

THET 

T) 

LAST 

AST  =  .13) 
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APPENDIX  F 


TEMPFRATURE    AND 
DATA 


RESISTIVITY    VALUES 
SET    A 


HEATING    CONDITIONS 


THE    NUMBER    OF    DATA    POINTS     IS  5 

THE    MAXIMUM    DEGREE    OF     FIT     IS  1 

THE    Y-AXIS    SCALE     IS  0.10    IN    DEGREES    PER    INCH 

THE    X-AXIS     SCALE     IS  1.00     IN    OHM-CM    PER     INCH 

THE    INITIAL    TEMPERATURE     (TO)     WAS    21.00IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NONDIMENSIONAL    VALUES 
THETA  RHO-PRIME 


27.5 

2753.6 

0.309524 

l.OOOOOQ 

31.3 

2508.5 

0.490476 

0.910969 

37.5 

2508.5 

0.7857  14 

U .9 10969 

44. C 

2014.9 

I. 095238 

0.731733 

51. 4 

1800.1 

1.447619 

0.653726 

COEFFICIENTS    OF    THE    POWER    SERIES    EXPANSION 
Y(X)=B(1)+6(2)*X+B(3 )*X**2. .. 

R(     1)=    1.09146D    00       B(     2) =-3. 02738D- 01 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 


FRRB(     1)=9. 2330-02       ERR6( 
SUM     SO    DEV    =    1.023D-03 
CHISO=       5.993D-03 


I 

X(I) 

F2(  I) 

1 

0.310 

1-00000000 

00 

2 

0.490 

9.1098925D- 

-01 

3 

0.786 

9. 1G98925D- 

-01 

4 

1.095 

7.3173300D- 

-01 

5 

1.448 

6-  5372603D- 

-01 

2)=1.008D-01 

F-RATIC  =  9.028D 
DEG  OF  FREEDOM  =   3. 


00 


Y(I) 

9.9775617D-01 
9.4297694D-01 
3.53  5  97  04D-01 
7.5989231D-01 
6.5321307C-01 


DELYU  ) 
-2.2418303D-03 
3.1987691D-02 

-5. 73922080-02 

2.8159305D-02 

-5 .12956950-04 
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TEMPERATURE  AND 
DATA 


RESISTIVITY 
SET  A 


VALUES 


COOLING    CONDITIONS 


THE  NUMBER    OE    DATA    POINTS     IS  9 

THE  MAXIMUM    DEGREE    OF     FIT     IS  1 

THE  Y-AXIS     SCALE     IS  0.10     IN    DEGREES    PER     INCH 

THE  X-AXTS    SCALE     IS  1.00     IN    OHM-CM    PER     INCH 

THE  INITIAL     TEMPERATURE     (TO)     WAS    21.00IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NONDIMENSIONAL    VALUES 
THETA  RHC-PRIME 


25.2 

2870.5 

0.200000 

1.000000 

29!. 9 

2687.9 

0.423810 

0.936367 

31.2 

2645.9 

0.485714 

0.921756 

32.7 

2546.3 

0.557143 

0.887058 

33!.  6 

2508.5 

0.600000 

0.873890 

34.6 

2453.8 

0.647619 

0.854834 

35.9 

2418.7 

0.709524 

0.842606 

40.8 

2242.2 

0.942857 

0.781113 

44.7 

2102.7 

1.128571 

0.732520 

COEFFICIENTS    OF     THE     POWER    SERIES    EXPANSION 
Y(X)=B(1)+B(2) ►X+B(3)? X**2... 

B(     1)=     1.05463D    00       B(     2) =-2. 9 1 728D-01 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 

FRRB(     1)=1. 9280-02       ERRB(     2)=2.313D-02 
SUM    SO    DEV    =    4.1700-05  F-RATIO    =     1.071D    02 

CHISQ=       4.3520-04  DEG    OF    FREEDOM    =       7. 


I 


X(I) 


F?  (  I) 


Y(  I) 


DEL  Yd  ) 


0.200       1.0000U00D    00       9. 96 279 8 40-01    -3 . 72 01594D-03 


2 

0.424 

9. 36387390-01 

3 

0.486 

9.217  55  7  9D- 

-01 

4 

0.5  57 

8.67C5800D- 

-01 

5 

0.600 

8.73889570- 

-Oi 

6 

0.648 

6.  5483365D- 

-01 

7 

0.710 

8*42605820- 

-01 

v. 309S636D-G1  -5 

9. 129290 1D-01  -8 

8.920913ID-01  5 

8.7958868D-01  5 

6 ,65696o6 D-01  1 

8.47637530-01  5 


3990293D-03 
82677790-03 
J033D-03 
69911680-03 
08632260-02 
0317150D-03 


0.943       7.8U1327D-01       7  .  79  56769D-01    -1 .  550  581  ID- 03 
1.129       7.3252047D-01       7. 253896 5D-01    -7 . 13081 33D-03 
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TEMPFRATJRE  AND 
DATA 


RESISTIVITY 
SET  B 


VALUES 


HEATING    CONDITIONS 


THE    NUMBER    OF    DATA    POINTS     IS  5 

THE    MAXIMUM    DEGREE     OF    FIT     IS  1 

THE    Y-AXIS    SCALE     IS  0.10    IN    DEGREES    PER    INCH 

THE    X-AXIS     SCALE     IS  1.00     IN    OHM-CM    PER     INCH 

THE    INITIAL    TE-MPERATURE     (TO)     WAS    23.00IN    DEG    C 


OUMENSl ONAL    VALUES 
TEMPER.ATURE       RESISTIVITY 
(DEGR&ES    C)  (OHM-CM) 


3723.4 
3422.1 

2997.6 
27  5  3.6 

2500.4 


NGNOIMENSIONAL    VALUES 
THETA  RHO-PRIME 


26 

.2 

32< 

.2 

40 

.4 

46 

.0 

52 

.0 

C0EFF1CI 

ENTS 

Y(  X) 

0.139130 
0.400000 
0.75  65  2  2 
1.000000 
1.260870 


OF    THE    POWER    SERIES    EXPANSION 
=B(1)+B(2)  <X+B(3)-Xw2.  .  . 

B(    1)=    1.03617D    00       B(     2 ) =-2. 93393D-0 1 

ESTIMATES    OF     ERROR    FOR    THE    COEFFICIENTS 


FRRB(     1 )=1 .533D-02        ERRB( 
SUtt     SO    ~,cV    =     3.4243-05 
CHISO=       2.218D-04 


F2(  I  ) 
1.00000000  00 
9.  1907934D-01 
8.O507063D-01 
7.3953913D-01 
6*7368534D-01 


2)=1.875D-02 

F-RATIO  =  2.449D 

DEG  GF  FREEDOM  =   3. 


1.000000 
0.919079 
0.805071 
0. 739539 
0.673685 


02 


I 

X(  I) 

1 

0.139 

2 

0.400 

3 

0.757 

4 

1.000 

5 

1.261 

Y(I)  DELY(I) 

9.9534654D-01  -4. 

9.1880931D-01  -2. 

8.14203440-01  9. 

7.4277369D-01  3. 

6.6623647D-01  -7. 


65346390-03 
7002533D-04 
1373018D-03 
2345ol6D-03 
4488742D-03 
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TEMPFRATURE    AND 
DATA 


RESISTIVITY    VALUES 
SET    B 


COOLING    CONDITIONS 


THE  NUMBER    OF    DATA    POINTS  IS 

THE  MAXIMUM    DEGREE    OF     FIT  IS 

THE  Y-AXIS     SCALE     IS  0.10     IN 

THE  X-AXIS    SCALF     IS  1.00     IN 


19 
1 

DEGREES    PER    INCH 
OHM-CM    PER     INCH 


THF     INITIAL    TEMPERATURE     (TO)     WAS    23.00IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NONDIMENS IONAL    VALUES 
THETA  RHO-PRIME 


24.1 

3872.4 

0.047826 

1.000000 

24.7 

3850.4 

0.073913 

0.994319 

25.7 

3764.8 

0.117391 

0.972214 

26.5 

3743.9 

0.152174 

0.966816 

2a. 3 

3662.9 

0.230435 

0.945899 

20.4 

3585.3 

0.278261 

0.925860 

30.0 

3547.7 

0.304348 

0.916150 

31.0 

3492. S 

0  .547326 

0.901973 

32.5 

3394.7 

0.413043 

0.87bo40 

35.7 

3226.3 

0.552174 

0.833153 

3  8.0 

3136.6 

0.652174 

0.809989 

40.3 

3079.5 

0.752174 

0.795243 

42.0 

2971.3 

0.826087 

0. 767302 

46.2 

2822.5 

1  .008596 

0.728876 

4C.0 

2  7  5  3  •  o 

1.086957 

0.711084 

51.4 

2546.2 

1.234783 

0.657i>25 

53.5 

2527.2 

i. 326087 

0.652619 

56.3 

2453.8 

1.447826 

0.633664 

5  9.6 

2384.6 

i  .591304 

0.615794 

COEFFICIENTS    OF    THE    POWER    SERIES    EXPANSION 
Y( X)=B(1)+B(2)*X+B(3)*X**2. . . 

R(    1)=    9. 968530-01       8(     2 ) =-2 . 59982D-01 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 

FRRB(     11=2.4040-02       ERRB(    2)=2. 9470-02 
SUM    SO    DEV    =     1.8400-04  F-RATIO    =    7. 7840    01 

CHISO=       4. 7870-03  DEG    OF    FREEDOM    =     17. 


I 
1 
2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1  3 

14 

1  5 

16 

17 

18 

19 


X(I) 

0.048 

0.074 

0-117 

0.152 

0.230 

0.27  3 

0.304 

0.343 

0.413 

0.552 

0.652 

0.75  2 

0.826 

1.009 

1.087 

1.235 

1.326 

1.443 

1.591 


F2(I  ) 
1  .00000 
9. 94318 
9.72213 
9.66816 
9.45899 
9-25859 
9.16150 
9. 01972 
8.!76639 
8.33152 
8.:09988 
7.!95243 
7.67301 
7<.28376 
7. 11083 
6.57525 
6..  52618 
6.33663 


6. 15793 


30D  00 
770-01 
61D-01 
4-4D-01 
18D-01 
93D-01 
19D-01 
940-01 
3 1 0-0 1 
570-01 
64D-01 
26D-01 
93D-01 
15D-01 
570-01 
05D-01 
53D-01 
680-01 
82D-01 


Y 

9.34 
9.77 
9.66 
9.57 
9.36 
9.24 
9.17 
9.06 
8.69 
8.53 
8.27 
8.01 
7.82 
7.34 
7.14 
6.75 
6.52 
6.20 
5.33 


(I) 

4186 
6364 
3329 
2900 
9436 
5097 
7275 
4240 
4686 
29  7  2 
2990 
3008 
0848 
6098 
2634 
8313 
0936 
4438 
1421 


2D-01 
9  0-01 
2C-01 
7D-01 
5D-01 
3D-01 
9  0-01 
30-01 
60-01 
8D-01 
8D-01 
SD^Ol 
2D-01 
6D-01 
4  0-01 
2D-01 
4D-01 
6D-01 
OC-01 


DELY 
1.558 
I.066 
5.  880 
9.526 
3.955 
1.350 
.577 
.451 
.282 
.014 
.  731 
.057 
.478 
.733 
.179 
.  630 
.246 
.322 
.265 


(I) 

13  76D-02 
223+D-02 
6934D-03 
3750D-03 
52990-03 
1984D-03 
H039D-03 
09410-03 
66750-02 
4712D-02 
0444D-G2 
6240D-03 
2693D-02 
70620-03 
87410-03 
62730-02 
5423D-04 
0020D-02 
17270-02 
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TEMPERATURE    AND    RESISTIVITY    VALUES 
DATA    SET    C 


HEATING    CONDITIONS 


THE  NUMBER    OF    DATA    POINTS    IS  5 

THE  MAXIMUM    DEGREE    OF    FIT     IS  1 

THE  Y-AXIS     SCALE     IS  0.10     IN    DEGREES     PER     INCH 

THE  X-AXIS     SCALE     IS  1.00     IN    OHM-CM    PER     INCH 

THE  INITIAL    TEMPERATURE    (TO)    WAS    22.50IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NONDIMENSIONAL    VALUES 
THETA  RHO-PRIME 


24.8 

3428.4 

0.102222 

1.000000 

29.5 

3110.6 

0.3111) .  1 

0.907304 

35  .4 

2799.1 

0.573333 

0.81644;> 

41.3 

2544.3 

0.835556 

0.742125 

47.8 

2299.9 

1.124444 

0.670838 

COEFFICIENTS    OF    THE    POWER    SERIES    EXPANSION 
Y( X  I sfl { 1 ) +B ( 2 )*X  +  B ( 3 )*X**2. . - 

R(     1J=     1.01512D    00       8(     2)=-3. 18632D-01 

ESTIMATES    OF    ERROR    FCR    THE    COEFFICIENTS 


FRRR(     1)=3. 2560-02       ERRB( 
SUM    SC    OEV    =     1. 7550-04 
CHISCJ=        1.0590-03 

I  X(I>  F2(I) 

I  X(I)  F2(  I  ) 

1  G.102  1 .OOOOOOOD  00 

2  0.311  9. 07303700-01 

3  0.573  8. 16444990-01 

4  0.836  7.4212461D-01 

5  1.124  6.7083771D-01 


2)=4. 701D-02 

F-RATIO  = 
DEG  OF  FREEDOM  = 


4.595D 
3. 


01 


Y(  I) 

Y<n 

9. 82551580-01    -1 
9.1599281D-01       8. 
8.3244032D-01       1. 
7.486tt783C-01       6. 
6. 5o83847D-01    -1 . 


DELYU  ) 
DELYU  ) 
7448425D-02 

68911200-03 
5995330D-02 
7632205D-03 
39992370-02 
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TEMPERATURE    AND 
DATA 


RESISTIVITY    VALUES 
SET    C 


COOLING    CONDITIONS 

THE    NUMBER    OF    DATA    POINTS     IS  21 

THE    MAXIMUM    DEGREE    OF     FIT    IS  1 

THE    Y-AXIS    SCALE     IS  0.10    IN    DEGREES    PER    INCH 

THE    X-AXIS    SCALE     IS  1.00     IN    OHM-CM    PER    INCH 

THE    INITIAL    TEMPERATURE    (TO)     WAS    22.50IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPFRATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NONDIMENS IONAL    VALUES 
THETA  RHO-PRIME 


23.5 

3463.8 

0.044444 

1.000000 

25.? 

3428.4 

0.120000 

0.989780 

25.8 

3359.8 

0.146667 

0.969975 

26.? 

3359.8 

0. 164444 

0.969975 

27.1 

3326.5 

0.204444 

0.960361 

28.4 

3261.8 

0.262222 

0.941683 

30.0 

3199.6 

0.333333 

0.923725 

31.0 

3169.4 

0.377778 

0.915007 

31.6 

3082.1 

0.404444 

0.689803 

33.3 

3026.4 

0.480000 

0.873723 

35.3 

2921.0 

0.568889 

0.843293 

37.2 

2822.7 

0.653333 

0.814914 

39.2 

2730.7 

0.742222 

0.788354 

40.8 

2  644.6 

0.813333 

0.763497 

46.7 

2415.9 

1.075556 

0.697471 

49.5 

2315.8 

1.200000 

0.668572 

52  .  0 

2253.5 

1.311111 

0.650586 

53.9 

2194.5 

1.395556 

0.633553 

55.6 

2  152  „  2 

1.471111 

0.621341 

57.5 

2085.2 

1.555556 

0.6  01998 

61.0 

1986.3 

1.711111 

0.573445 

COEFFICIENTS    OF    THE     POWER    SERIES    EXPANSION 
Y  (  X  )  =  B  (  1 )  +  B  (  2  J  -■■=  X+  B  (  3  )  *  X**  2 .  .  . 

R(     1)=    I.OO0O8D    00       B(     2) =-2.634730-01 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 
FRRR(      1)=2. 3390-02        ERR6(     2)=2.634D-02 
SUM    SO    DEV    =     1.7330-04  F-RATIO    =     1.0390    02 

CHISO=       5.0500-03  DEG    OF    FREEDOM    =     19. 


I 

xm 

1 

0.044 

2 

0.120 

3 

0.147 

4 

0.  164 

5 

0.204 

6 

0.262 

7 

0.333 

8 

0.378 

9 

0.404 

1  0 

0.480 

11 

0.569 

12 

0.653 

13 

0.742 

14 

0.813 

15 

1.076 

16 

1.200 

17 

1.311 

1  8 

1.396 

19 

1.47  1 

20 

1.55  6 

21 

1.711 

F2  (  I  ) 

1.00000000  00 
9. 89780010-01 
9.:6997517D-01 
9.'699751  70-01 
9. 60361450-01 
9,41682550-01 
9,23725390-01 
9. 15006640-01 
8.8980311D-01 
8.7372250D-01 
8u43293490-01 
8. 1491426D-01 
7.88353830-01 
7.63496740-01 
6. 97470990-01 
6.:6857209D-01 
6.50586060-01 
6.3355275D-01 
6-21340720-01 
6.01997310-01 
5. 73445350-01 


Yd) 

9.94148460-01 
9.7386386D-01 
9.6670459C-01 
9. 6 193 174  0-01 
9.5119283D-01 
9.3568107D-01 
9.1658968D-01 
9.04657560-01 
8.97498290-01 
8. 77213690-01 
8.5334945D-01 
8.30678420-01 
8.06814180-01 
7. 87722790-01 
7.17323280-01 
6.83913340-01 
6.5408305D-01 
6.31412C2D-01 
6. 11127410-01 
5.8845639D-G1 
5. 46693970-01 


DELY(I) 

-5 .85153760-03 

-i .59161520-02 

-3.27058510-03 

-8 .04343310-03 

-9.  1686220D-03 

-6.00147030-03 

-7.13570530-03 

-1. 03490770-02 

7. 69518480-03 

3 .4911661D-03 

1.0055955D-02 

1 .57641640-02 

1. 8460349D-02 

2.42260500-02 

1.98522950-02 

1. 53412560-02 

3.49698360-03 

-2.14072C2D-03 

-1.0213311D-02 

-1 .35414  200-0  2 

-2.67513830-02 
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TEMPERATURE    AND 
DATA 


RESISTIVITY    VALUES 
SET    D 


HEATING    CONDITIONS 


THE  NUMBER    OF     DATA    POINTS     IS  5 

THE  MAXIMUM    DEGREE    OF     FIT     IS  1 

THE  Y-AXIS     SCALE     IS  0.10    IN    DEGREES    PER    INCH 

THE  X-AXIS     SCALE     IS  1.00     IN    OHM-CM    PER    INCH 

THE  INITIAL     TEMPERATURE     (TO)     WAS    19.00IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(OEGREES    C)  (OHM-CM) 


24.7 

32.2 
40.3 
46.5 
54.8 


4271.7 
3551.6 
3182.6 
2&10.8 
2  49  6.5 


NONDIMENS IONAL  VALUES 
THETA        RHO-PRIME 


0.300000 
0.694737 
1.121053 
1.447368 
1.884211 


1.000000 
0.831425 
0.745043 
0.658005 

0. 584428 


COEFFICIENTS    OF    THE    POWER    SERIES 
Y(X)  =8(1  ) +6(2  )-:X  +  b(  3  )--X- 


EXPANSION 
■2  . .  . 


R(     1)=     1.04268D    00       b(     2 ) =-2 . 55 9960-0 1 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 

FRRB(     l)=7.l67D-02       ERRB(     2)^5.8620-02 
SUM    SO    DEV    =     6.3490-04  F-RATIC    =    1.9070    CI 

CHIS0=       3.9800-03  DEG    OF     FREEDOM    =        3. 


I 

X(  I) 

F2(I) 

1 

0.30  3 

I. 

.00000000 

00 

2 

0.695 

8, 

,31425430- 

-01 

3 

1.121 

7. 

,145042960- 

-01 

4 

1.447 

6. 

.580050  1D- 

-01 

5 

1.884 

5. 

.&442775D- 

-01 

Y(I) 

9.6588252D-01 
8. 6463137 C-01 

7.5569614D-01 
6.7216052D-01 
5. 60330590-01 


DELYCI ) 

-3. 41174820-02 
3.34059^50-02 
1. 0653179D-02 
1.41555140-02 

-2.4097156D-02 
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TEMPERATJRE    AND 
DATA 


RESISTIVITY    VALUES 
SET     D 


COOLING  CCMDITIONS 


THE  NUMBER    OF    DATA    POINTS  IS 

THE  MAXIMUM    DEGREE    OF    FIT  IS 

THE  Y-AXIS     SCALE    IS  0.10     IN 

THE  X-AXIS     SCALE     IS  1.00    IN 


19 
1 

DEGRFES    PER    INCH 
OHM-CM    PER    INCH 


THE     INITIAL    TEMPERATURE     (TO     WAS     19.00IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NONDIMENSIONAL    VALUES 
THETA  RHO-PRIME 


24.1 

4166.1 

0.268421 

1.000000 

24.0 

4166. 1 

0.263158 

1.000000 

25.5 

4065.7 

0.342105 

0.975901 

25.7 

4017.2 

0.352632 

0.964259 

28.2 

3923.7 

0.484211 

0.941816 

29.9 

3791.3 

0.573684 

0.910036 

31.0 

3749.1 

0.631579 

0.899906 

33.5 

3551.6 

0.763158 

0.852d00 

35.7 

3407.8 

0.873947 

0.817983 

36.3 

3442.7 

1.015789 

0.826360 

42.5 

3275.4 

1.236842 

0.786203 

4b.. 4 

3066.7 

1.442105 

0.736108 

47.6 

3066.7 

1.515789 

0.736108 

50.4 

2498.0 

1.652632 

0.599602 

52.0 

2461  .5 

1.73  6  34  2 

0.590840 

54.7 

2408.6 

1.878947 

0.578143 

55  .6 

2358.0 

1.936842 

0.565997 

55.3 

2189.2 

1.910526 

0.525479 

57.8 

2120.3 

2.042105 

0.508941 

COEFFICIENTS    3 F     THE    PCW=R    SERIES    EXPANSION 
Y(X)=B(1 )+6<2)*X+B(3)*X**2.. . 

R(     1>  =     1.07062D    00       B(     2) =-2. 640460-01 

ESTIMATES    GF    ERROR    FOR    THE    COEFFICIENTS 

FRRB(     11=5.7410-02       ERRB(     2)=4. 5320-02 
SUM    SO    DEV    =    7.1970-04  F-RATIO    =     3.3950    01 

CH1S0=       2.0390-02  DEG    OF    FREEDOM    =    17. 


I 

X(I) 

F2U) 

Y(  I) 

1 

DELYd  ) 

1 

0.268 

1. 

OOOOOOOD 

00 

9. 

99739700-01 

-2, 

.60300800- 

•04 

2 

0.263 

1. 

OOOOOOOD 

00 

1. 

00112940 

00 

1, 

. 12941360- 

-03 

3 

0.342 

9. 

■  75  9  0072  0- 

-01 

9. 

,30283700- 

-01 

4 

.3829747D- 

•03 

4 

0.353 

9. 

64259140- 

-01 

9. 

7750427D- 

-01 

1, 

, 3245129D- 

-02 

5 

0.464 

9. 

,41816090- 

-01 

9. 

,42761410- 

•01 

9, 

. 45320630^ 

■04 

6 

0.574 

9. 

1003576D- 

-01 

9. 

19136260- 

-Oi 

9 

.10049  75  0- 

•03 

7 

0.632 

8. 

,9990639D- 

-01 

Q 

,  03849<tOD-01 

3, 

.94301640- 

•03 

8 

0.763 

8- 

,5?49994D- 

-01 

8, 

,69106540- 

-01 

1 

. 66066G3D- 

•02 

9 

0.879 

8* 

179S325D- 

-01 

8. 

38532830- 

-01 

2 

.05495800- 

-02 

10 

1.016 

8. 

,:2636039D- 

-01 

8. 

,02400250- 

-01 

-2, 

.39601350- 

•02 

1  1 

1.237 

7- 

.662  02  92  D- 

-01 

7, 

,44032240- 

-01 

-4 

.21706790- 

-02 

12 

1.442 

7. 

,3610611D- 

-01 

6. 

,89833380- 

-01 

-4. 

.62747290- 

-02 

13 

1.516 

7. 

,36108110- 

-01 

6. 

,70377380- 

-01 

-6 

.573C731D- 

•02 

14 

1.653 

5„ 

99601550- 

-01 

6. 

342t480D- 

-01 

3 

.46432590- 

-02 

15 

1.73  7 

5i 

,  90840350- 

-01 

6. 

,  120093  70- 

-01 

2, 

. 11690200- 

-02 

16 

1.879 

5. 

.78142  6  3  0- 

-01 

5, 

,74487080- 

-01 

-3 

.65554310- 

■03 

1  7 

1.937 

5. 

,c5996980- 

-01 

5. 

5920023D- 

-01 

-6 

.79675000- 

-03 

18 

1.911 

5. 

,25479470- 

-01 

5. 

66148300- 

-01 

4, 

.06693330- 

■02 

19 

2.042 

5. 

,08941220- 

-01 

5. 

,31405940- 

-01 

2 

.2464721D- 

•02 
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TEMPERATURE    AND    RESISTIVITY    VALUES 
DATA    SET     E 


HEATING    CONDITIONS 


THE  NUMBER    OF    DATA    POINTS     IS  5 

THE  MAXIMUM    DEGREE    OF     FIT     IS  1 

THF  Y-AXIS     SCALE    IS  0.10    IN    DEGREES    PER    INCH 

THE  X-AXIS     SCALE     IS  1.00    IN    OHM-CM    PER     INCH 

THE  INITIAL    TEMPERATURE    (TO)     WAS    22. 50  IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


N0ND1MENSI0NAL    VALUES 
THETA  RHO-PRIME 


25.7 

3451.0 

0.1422  22 

1.000000 

35.7 

310  5.5 

0.586667 

0. 899884 

41.0 

2894.7 

0.622222 

0.8  36  600 

44.7 

2710.6 

0.986667 

0.785453 

49.5 

2529.6 

1.200000 

0. 733005 

COEFFICIENTS    OF    THE     POWER    SERIES    EXPANSION 
Y( X ) =  B ( 1  J  +B ( 2 ) *X  +  B ( 3 )*X**2 . . . 

R(     1)=     1.04182D    00       B(    '2)  =  -2.  54679D- 01 

.ESTIMATES    OF     ERROR    FOR    THE    COEFFICIENTS 

FRRB<     l)=1.6Q2D-02        ERRB(     2)=2.036D-02 
SUM     SO    DEV    =    3.2810-05  F-RATIC    =     1.5650    02 

CHISC=        1.S93D-04  DEG    OF    FREEDOM    =       3. 


I 

XI 

I) 

F?(  I  ) 

1 

0. 

.142 

UOOOOOOOD 

00 

2 

0. 

.587 

8.9986409D- 

-01 

3 

0. 

.822 

6.3680035D- 

-01 

4 

0. 

.987 

7.6545349D- 

-01 

5 

1. 

.200 

7.3300493D- 

-01 

Y(I) 

1.0055945D  00 
8.92403t>7D-01 
8.3241250D-01 
7. 90531680-01 
7.3620026D-01 


DELY(I) 

5.5945479D-03 

-7.4804221D-03 

-6.3878439D-03 

5.0783671D-03 

3. 19533100-03 
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TEMPERATURE    AND 
DATA 


RESISTIVITY    VALUES 
SET    E 


COOLING    CONDITIONS 


THE    NUMBER    OF    OATA    POINTS     IS  9. 

THE    MAXIMUM    DEGREE    OF     FIT    IS  1 

THE    Y-AXIS     SCALE     IS  0.10    IN    DEGREES    PER     INCH 

THE    X-AXIS    SCALF     IS  1.00     IN    OHM-CM    PER     INCH 

THE     INITIAL     TEMPERATURE    (TO)     WAS    22.50IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


3796.6 
3674.0 
3559.0 
3522.3 
3416.4 
3382.6 
3316.8 
3284.9 
2627.1 


NCNOIMENSIONAL    VALUES 
THET*.  RHO-PRIME 


24, 

.3 

26 

.4 

27. 

.5 

29, 

.5 

30, 

.5 

32, 

.0 

33 

.0 

33, 

.3 

63, 

►  5 

COEFFICIENTS 

Y(X) 

0.080000 
0.173333 
0.222222 
0.311111 
0.355556 
0.422222 
0.466667 
0.480000 
1.822222 


1.000000 
0.967708 
0.937418 
0.927751 
0.89985S 
0.890955 
0. 873624 
0.865222 
0.691961 


OF    THE     POWER    SERIES    EXPANSION 
;  6  (  1 )  +B  (  2  )  '  X+  B  (  3  )  •■  X**2.  .  . 

R(     1)=    9.73366D-G1       B(     2) =-1. 62 877D-01 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 


FRRBi     l)  =  3.445D-02       ERRBi 
SUM    SO    DEV    =    4. 7060-04 
CHISC=       4.65.7C-03 


F?(  I  ) 

1.00000000  00 
9. 67707950-01 
9.3741 76 9D-01 
9..2775115D-01 
8.99857770-01 
8.9G9  55  07D-ul 
8. 7362377D-01 
8.6522151D-01 
6.:9196123D-01 


I 

xm 

1 

0.080 

2 

0.173 

3 

0.222 

4 

0.311 

5 

0.3  56 

6 

0.422 

7 

0.467 

8 

0.480 

9 

1.822 

2) =5. 010D-02 

F-PATIO 
DEG  OF  FREEDOM 

Yd) 

9.6033601D-01 
9.4513418C-01 
9.3  7171320-01 

9.226SS390-01 
9.  1545442^-01 
9.0459597D-01 
8.9735700D-01 
8.9518531D-01 
6.76568530-01 


1.057D 
7. 


01 


DELYU  ) 
•3  .96639890-02 
2.2573768D-02 
2 .46371 19D-04 
•5. 05776030-03 
1.55966520-02 
1.36  4090  40-0  2 
2.3733234D-02 
2.99637990-02 
1.5392701C-02 
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TEMPERATURE    AND    RESISTIVITY    VALUES 
DATA    SET     F 


HEATING    CONDITIONS 


THE  NUMBER    OF    DATA    POINTS    IS  5 

THE  MAXIMUM    DEGREE    OF     FIT     IS  1 

THE  Y-AXIS    SCALE     IS  0.10     IN    DEGREES    PER    INCH 

THE  X-AXIS     SCALE     IS  1.00     IN    OHM-CM    PER     INCH 

THE  INITIAL    TEMPERATURE    (TO)    WAS    22.5DIN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE        RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NGNDIMENSIONAL    VALUES 
THETA  RHO-PRIME 


24.2 

3298.8 

0.075556 

1.000000 

30.0 

3009. 1 

0.  333333 

0.912180 

33-. 0 

2559.3 

0.686389 

0.775828 

44.8 

2381.3 

0.991111 

0.721669 

49.0 

2197.8 

1.177776 

0.666242 

COEFFICIENTS    DF    THE    PGWER    SERIES    EXPANSIGN 
Y(X)=8(1 )+B(2)*X+B(3)*X**2.. . 

B(     1)=     1.011910    00      B(     2) =-3.010450-01 

ESTIMATES  OF  ERROR  FOR  THE  COEFFICIENTS 

FRRBt  J)=3. 6070-02   ERRB(  2>=4. 6870-02 

SUM  SO  DEV  =  2.1790-04  F-RATIO  =  4.1250  01 

DEG  OF  FREEDOM  =   3. 


CHIS0=   1 

.3600-03 

I    X(I) 

F2(I) 

1    0.076 

1..00000000 

00 

2    0.333 

9.1218  019D- 

-01 

3    0.689 

7.75827570- 

-OL 

4    0.991 

7.2186856D- 

•01 

5    1.178 

6 .66242270- 

-0  1 

Yd ) 

9.8916102D-01 

9. 11558220-01  -6 

6.0^5196  8  0-01  2 

7. 13537300-01  -8 

6.5734217D-01  -8 


DELYU  ) 
1. 08389630-02 
21964870-04 
86923100-02 
33126160-03 
90010070-03 
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TEMPERATURE  AND 
DATA 


RESISTIVITY 
SET  F 


VALUES 


CODLING    CONDITIONS 


THE  NUMBER    OF    DATA    POINTS     IS  12 

THE  MAXIMUM    DEGREE    OF     FIT     IS  1 

THE  Y-AXIS     SCALE     IS  0.10    IN    DEGREES    PER    INCH 

THE  X-AXIS     SCALF     IS  1.00     IN    OHM-CM    PER     INCH 

THE  INITIAL    TEMPERATURE     (TO)     WAS    22.50IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NONDIMENS IONAL    VALUES 
THETA  RHO-PRIME 


24.0 

3630.9 

0.066667 

1.000000 

25.0 

3581.6 

0.111111 

0.986422 

25.5 

3574. L 

0.133333 

0.984356 

27.9 

3465.7 

0.240000 

0.954502 

29.8 

3330.9 

0.324444 

0.917376 

33.5 

3206.2 

0.488839 

0.883032 

38.2 

2906.9 

0.697778 

0.800600 

43.0 

2743.9 

0.911111 

0.755708 

51.0 

2449.4 

1.266667 

0.674599 

55.0 

2348.6 

1.444444 

0.646837 

60.5 

2183.8 

1.688889 

0.601449 

65.5 

2077.7 

1.911111 

0.572227 

COEFFICIENTS    OF    THE     POWER    SERIES    EXPANSION 
Y(XI=B( 1)+B( 2)*X+B(3)*X#*2... 

B<     1)=    L.OO2O0D    00       B(     2) =-2. 42002D-01 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 

FRRBt     l)=3.175C-02       ERRB(     2)=3.181D-02 
SUM    SO    DEV    =    3. 3503-04  F-RATIO    =     5.786D     01 

CHISQ?       5.607D-03  DEG    OF    FREEDOM    =    10. 


I 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1  1 

12 


X(I) 

C.06  7 

0.111 

0.133 

0.240 

0.324 

0.489 

0.698 

0.911 

1.267 

1.44  4 

1.68  9 

1.911 


F2 
1..00 

9..  86 
9.:34 
9*  54 
9.17 
8-8-3 

8.:oo 
7.:55 
6.74 
6.46 
6-OL 
5-72 


(  I  ) 

0000 
4221 
3565 
5016 
3758 
0317 
6004 
7079 
5985 
8368 
4486 
2272 


OD  00 
0D-01 
OD-01 
40-01 
6D-01 
6D-01 
OD-01 
5D-01 
8D-01 
7D-01 
3D-01 
7D-01 


Y 
9.85 
9.75 

9.69 
9.^3 
9.23 
8.83 
8.33 
7.81 
6.95 
6.  52 
5.93 
5.39 


(I) 

86368D-01 
10301D-01 
73018D-01 
916580-01 
48082D-01 
68486D-01 
133230-<01 
506030-01 
4607  1D-01 
43605D-01 
28189D-01 
50356D-0I 


DELY 
.413 
.  131 
.  462 
.058 
.104 
.531 
.253 
.579 
.086 
.601 
.  166 
.2  72 


(I  ) 

63220-02 

4087D-02 

631o0-02 

50560-02 

96100-03 

00990-04 

2827D-02 

80840-02 

2127D-02 

1774D-03 

78660-03 

37080-02 
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TEMPFRATURE  AND 
DATA 


RESISTIVITY 
SET  G 


VALUES 


HEATING    CONDITIONS 


THE    NUMBER    OF     DATA    POINTS     IS  5 

THE    MAXIMUM    DEGREE    OF     FIT    IS  1 

THE    Y-AXIS    SCALE    IS  0.10     IN    DEGREES    PER     INCH 

THE    X-AXIS     SCALC     IS  1.00     IN    OHM-CM    PER     INCH 

THE    INITIAL    TEMPERATURE    (TO)     WAS    15.00IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


25.0 
30.3 
39.0 
47.3 
51.4 


3134.0 
2  822.8 
2439.1 
2188.5 
2031.9 


NONDIMENSIONAL    VALUES 
THETA  RHO-PRIME 


0.666667 
1.0200  00 
1.600000 
2.153333 
2.426667 


1.000000 
0.900702 
0.778271 
0.698309 
0.648341 


COEFFICIENTS    OF    THE    POWER    SERIES     EXPANSION 
Y(  X)  =  S(  1)+B( 2)*X+B(3)*X**2... 

R(     1)=    1.11070D    00       B(    2)=-1.94223D-01 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 

FRRB(     1)^4.8860-02       ERRB(     2)=2.862D-02 
SUM    SO    DEV    =     2.156D-04  F-RATIO    =     4.605D    01 

CHISO=       1.27oD-03  DEG    OF    FREEDOM    =       3. 


1 

0.667 

l.'OOOOOOOD 

00 

2 

1.020 

9«0070198D- 

-01 

3 

1.600 

7...782  705  8D- 

•01 

4 

2.153 

6.983088  7 D- 

-01 

5 

2.42  7 

6*48340730- 

-01 

9.8122036D-01 
9.1259474D-01 
7.9994515D-01 
6.9247485D-01 
6.3938711D-01 


-1. 8779641D-02 

1 .18927660-02 

2. 16745690-02 

-5. 83402310-03 

-8.9536710D-03 
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TEMPERATURE    AND    RESISTIVITY    VALUES 
DATA    SET    G 


COOLING    CONDITIONS 


THE  NUMBER    OF    DATA     POINTS    IS  10 

THE  MAXIMUM    DEGREE    OF    FIT     IS  1 

THE  Y-AXIS    SCALE    IS  G.10     IN    DEGREES    PER    INCH 

THE  X-AXIS     SCALF     IS  1.00    IN    OHM-CM    PER    INCH 

THE  INITIAL    TEMPERATURE     (TO)     WAS     15.00IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NDNDIMENSIONAL    VALUES 
THETA  RHO-PRIME 


23.0 

3416.5 

0.533333 

1.000000 

24.7 

3253.6 

0.646667 

0.952320 

28.0 

3105.5 

0.866667 

0.908971 

30.7 

2970.3 

1.046667 

0.869399 

35.5 

2732.3 

1.366667 

0.799737 

39.0 

2587.2 

1.600000 

0.757266 

45.5 

2371.3 

2.033333 

0.694073 

50.8 

2217.0 

2.386667 

0.648910 

56.3 

2069.0 

2.753333 

0.605591 

63-2 

1917.5 

3.213333 

0.561247 

COEFFICIENTS    OF    THE     POWER    SERIES    EXPANSION 
Y(X)=3( 1)+B( 2)*X+B(3)*X**2... 

R(     1)=     1.04624D    00       R(     2)  =-1. 620310-01 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 


FRRB(      l)=5.445D-02       ERRB( 
SUM    SO    DEV    =     5.236D-04 
CHI SQ=       7. 1220-03 


I 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


xm 

0.  53  3 
0.647 
0.867 
1.04  7 
1.367 
1.600 
2.033 
2.387 
2.753 
3.213 


F2(I  ) 

1-00000 
9*52319 
9.08971 
8. :693  9b 
7.99736 
7,57266 
6.94072 
6i.  48909 
6. 05590 
5. .'61246 


OOD  00 
6  3D- 01 
17D-01 
51D-01 
570-01 
21D-01 
6  80*01 
70D-01 
52D-01 
89D-01 


2)=2. 922D-02 

F-RATIO    =    3.0  74D 
DEG     OF    FREEDOM    =       8. 


01 


9 
9 

9. 
8, 


Y 
59 
41 
05 
76 
8.24 
7.86 
7.  16 
6.59 
6.00 
5.25 


(  I) 

8211 

4577 

8109 

6454 

7957 

9885 

7753 

5245 

1133 

5792 


6D-01 
10-01 
9  0-01 
90-01 
1U-01 
7  0-01 
3D- 01 
3D- 01 
2D-01 
6D-01 


DEL 

■4.0  1 

■1.  08 

•3.16 

7.24 

2.50 

2.97 

2.27 

1.06 

■5  .47 

■3.  56 


Yd  ) 

768210-02 
61910D-02 
01795D-C3 
69802D-03 
591330-02 
22365D-02 
02445D-02 
148220-02 
719970-03 
676350-02 
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TEMPERATURE     AND 
DATA 


RESISTIVITY    VALUES 
SET    H 


HEATING    CONDITIONS 


THE  NUMBER    OF    DATA    POINTS    IS 

THE  MAXIMUM    DEGREE    OF    FIT     IS 

THE  Y-AXIS     SCALE     IS  0.10     IN 

THE  X-AXIS     SCALE     IS  1.00     IN 

THE  INITIAL    TEMPERATURE     (TO)     WAS 


5 
1 

DEGREES    PER    INCH 
OHM-CM    PER    INCH 
7.50IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE       RESISTIVITY 
(DEGREES    U  (OHM-CM) 


NONDIMENSIONAL    VALUES 
THETA  RHO-PRIME 


24.7 

3062.9 

2.293333 

1.000000 

3  0 .  0 

2  744.0 

3.000000 

0.895883 

40.7 

2316.8 

4.42b6to7 

0.756407 

45.4 

2142  .7 

5.0533  33 

0.699566 

50.7 

1970.0 

5.760000 

0-643181 

R  (     1 )  = 


COEFFICIENTS    OF     TH?     POWER    SERIES    EXPANSION 
Y(X)=B( 1»+3(2)*X+B(3)*X**2. .. 

i94D    00       B(     2) =-1.015260-01 


ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 


FRRB(     1)=5.413D-G2       ERRB(     2>=1.258D-02 
SUM    SO    DEV    =     1.566D-04  F- RATIO    =     6.5120    01 

CHISO=       9.317D-04  DEG    OF    FREEDOM    =       3. 


X(I) 

F2(  I  ) 

1 

2.29  3 

1.00000000 

00 

2 

3.00  0 

3-9588299D- 

-01 

3 

4.427 

7-56407330- 

-01 

4 

5.053 

6.99565770- 

-01 

5 

5.760 

6.43181300- 

-01 

Y(I) 

9.3310763D-01 
9.  1136272D-01 
7.6o5192  1D-01 
7.0239636D-01 
6.3115145D-01 


DELYd  ) 

-1 .68923650-02 
1.5479733D-02 
1.0111887D-02 
3.3305936D-03 

-1.2029849D-02 
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TEMPERATURE    AND 
DATA 


RESISTIVITY 
SET    H 


VALUES 


COOLING    CONDITIONS 


THE  NUMBER    OF    DATA    POINTS     IS 

THE  MAXIMUM    DEGREE    OF     FIT     IS 

THE  Y-AXIS     SCALE     IS  0.10    IN 

THE  X-AXIS     SCALE     IS  1.00    IN 

THE  INITIAL    TEMPERATURE     (TO)     WAS 


10 

1 
DEGREES    PER    INCH 
OHM-CM    PER     INCH 
7.50IN    DEG    C 


DIMENSIONAL    VALUES 
TEMPERATURE        RESISTIVITY 
(DEGREES    C)  (OHM-CM) 


NONDIMENSIONAL    VALUES 
THETA  RHO-PRIME 


24.0 

3267.4 

2.200000 

1 .000000 

24.5 

3267.4 

2.266667 

1.000000 

2  5.0 

3267.4 

2.333333 

1.000000 

26.5 

3147.3 

2.533333 

0.963243 

3L.4 

2901.0 

3.186667 

0. 887862 

37.2 

2658.7 

3.960000 

0.813705 

44.5 

2364.3 

4.933333 

0.723756 

49.5 

2226.4 

5.600000 

0.681398 

65.5 

1842.6 

7.733333 

0.563935 

70.5 

1757.4 

8.400000 

0.537859 

COEFFICIENTS    OF    TH^     POWER    SERIES    EXPANSION 
Y (  X )  =  3 ( 1 ) +3 ( 2 ) -X  +  3 ( 3  ) *X**2 . . . 

R(     1)=     1.15284D    00      B(     2 ) =- 7. 77955D- 02 

ESTIMATES    OF    ERROR    FOR    THE    COEFFICIENTS 

FRRB(     1)  =  7. 0090-02       ERRB(     2)=1.450D-02 
SUM    SO    DtV    =    7-9990-04  F-RATIC    =    2.860D    01 

CHISO=        1.C99D-02  '        DEG    OF    FREEDOM    =       8. 


I 

xm 

FPU  ) 

Y(  I) 

DELYU  ) 

1 

2.200 

1. 

0  000  00  00 

00 

9. 

816872  7D- 

-01 

-1.8312735D- 

-02 

2 

2.2o7 

1  . 

,00000000 

00 

9. 

,7b50090D- 

-01 

r-2.3499101r> 

-02 

3 

2.333 

U 

oooocooo 

00 

9, 

,  7  131453D- 

-01 

-2 .8  68  5467  0- 

-02 

4 

2.53  3 

9. 

,63242950- 

-01 

9. 

55755440- 

-01 

-7.46751030- 

-03 

5 

3.187 

6.1 

,67361910- 

-01 

9. 

,04929050- 

-01 

1.70671400- 

•02 

6 

3.960 

8. 

1370509D- 

-01 

8. 

4476720D- 

-01 

3.10621160- 

-02 

7 

4.933 

7. 

,23755390- 

-01 

7. 

6904626D- 

-01 

4. 5  2  903680- 

-02 

8 

5.600 

6. 

'8139805D- 

-01 

7. 

,17182600- 

-01 

3.57845460- 

-02 

9 

7.733 

5i. 

63934630- 

-01 

5. 

51218390- 

-01 

-1.2715739D- 

-02 

0 

8.40  0 

5. 

,37858650- 

-01 

4. 

99355230- 

-01 

-3. 85036200- 

-02 

106 


SAMPLE  DATA 
SPHERICAL  PROBE 

DATA  SET  I 

T  I  TAU  R  R/RO  R  R/RO  R  R/RO 

(SFC)        (AMPS)  (CM)  (CM)  (CM) 

R-27  ROCHROME  R-37 

30.0  C.050  0.78  1.05  1.62 

60.0  0.050  1.56  1.20  1.85            MOT                                  NOT 

120.0  0.050  3.12  1.45  2.23 

240.0  0.050  6.25  1.85  2.85            USED                               USED 

480.0  0.050  12.50  2.00  3.08 

960.0  0.050  25.00  2.30  3.54 

THFTA/GO  0.1260 

RESISTIVITY    (OHM-CM)  2420.0 

ISOTHERM    (DEGREE    C)  27.9 

THF     IMTIAL    TEMPERATURE    WAS     19.2     DEGREES    C. 

THP    PROSE    RADIUS    (ACTUAL)     WAS    0.241    CM. 

THF     PRCBE    RADIUS     IN    THE    PICTURES    WAS     0.650    CM. 


SPHERICAL    PROBE 


DATA     SET 

J 

T 
(SFC) 

I 

(AMPS  ) 

TAU 

R            R/RO 
(CM) 

R            R/RO 
(CM) 

R            R/RO 
(CM) 

30.0 
60.0 
120.0 
240.0 
480.0 
960.0 

0.060 

0.060 
0.060 
0.060 
0.060 
0.060 

0.78 
1.56 
3.12 
6.25 
12.5 
25.00 

R-27 
1.10       2.90 
1.40       3.68 
1.65       4.34 
2.15       5.66 
2.65       6.97 
3.30       8.68 

ROCHROME 
0.90       2.37 
1.10       2.90 
1.40       3.68 
1.70       4.47 
2.05       5.40 
2.55       6.71 

R-37 
0.0         0.0 
0.90       2.37 
1.00       2.63 
1.25       3.29 
1.45       3.82 
1.70      4.47 

THFTA/GO 

RESISTIVITY     (OHM-CM) 
ISOTHERM    (DEGREE    C) 

0.0200 
3605.0 
27.9 

0.0282 
3520.0 
29.9 

0.0572 
3190.0 
37.0 

THF     INITIAL    TEMPERATURE     WAS    23.0    DEGREES    C. 

THE     PROBE    RADIUS'    (ACTUAL)    WAS    0.241    CM. 

THF     PRCBE    RADIUS    IN    THE     PICTURES     WAS    0.38    CM 
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SAMPLE  DATA 
CYLINDRICAL  PROBE  (L/D=4) 

DATA  SET  K 

T  I  TAU  R  R/RO  R  R/RO  R  R/RO 

(SFC)        (AMPS)  (CM)  (CM)  (CM) 

R-27  ROCHROME  R-37 

30.0  0.030  3.11  0.70  3.68  0.0  0.0               NOT 

60.0  0.030  6.23  0.90  4.74  0.50  2.63 

120.0  0.030  12.45  1.15  6.05  0.80  4.21                USED 

240.0  0.030  24.91  1.42  7.47  1.10  5.79 

THFTA/GG  0.367  0.501 

RESISTIVITY     (OHM-CM)  3455.7  3324.9 

ISOTHERM    (DEGREE    C)  27.9  29.9 

THF     INITIAL    TEMPERATJRE    WAS    22.5  DEGREES    C. 
THF     PROBE    RADIUS     (ACTUAL)     WAS    0.121    CM. 

THE     PROBE    RADIUS    IN    THE    PICTURES  WAS    0.19    CM. 
THF    PROBE    LENGTH    WAS    0.96    CM. 


CYLINDRICAL    PROBE     (L/D=10) 

DATA    SET  L 

T                     I               TAU            R            R/RO  R             R/RO             R            R/RO 

(SFC)        (AMPS)                         (CM)  (CM)                            (CM) 

R-27  ROCHROME                     R-37 

60.0       0.090          6.23       0.80       4.00  0.5C       2.50                NOT 

120.0       0.090       12.45       1.20       6.00  0.90       4.50                   USED 

240.0       0.090       24.91       1.55       7.75  1.30       6.50 

THFTA/GO                                               0.-357  0.452 

RESISTIVITY    ICHM-CM)              3465.7  3324.9 

ISOTHERM    (DEGREE    f . )                      27.9  29.9 

THC     INITIAL    TEMPERATURE     WAS    20.4  DEGREES    C. 
THF     PROBE    RADIUS     (ACTJAL)     WAS    0.121    CM. 

THF     PROBE    RADIUS    IN    THE    PICTURES  WAS    0.200    CM. 
The    PROBE    LENGTH    WAS    2.42    CM. 


CYLINDRICAL    PROBE    (L/D=20) 

DATA    SET  M 

T                      I                 TAU             R             R/RO  R             R/RO              R             R/RO 

(SFC)       (AMPS)                        (CM)  (CM)                           (CM) 

R-27  ROCHROME                     R-37 

30.0       0.150         3.11      0.70      4.31  0.48      2.95          CO         0.0 

60.0       0.150          6.23       1.20       7.39  0.80       4.92          0.0          0.0 

120.0       C.150       12.45       1.62      9.97  1.15       7.08          0.45       2.77 

240.0       0.150       24.91       2.20    13.54  1.62       9.97          0.70       4.31 

THFTA/GO                                             0.356  0.489                          NOT 

RESISTIVITY     (GhM-CM)           3551.7  3480.0                           USED 

ISOTHFRM    (DEGREE    C)                   27.9  29.9 

THF     INITIAL     TEMPERATURE     WAS    22.5  DE3REES    C. 

THF     PRQBF    RADIUS     (ACTUAL)     WAS    0.121    CM. 

THF     PROSE    RADIUS    IN    THE    PICTURES  WAS    0.163    CM. 

THF    PRGBE    LENGTH    WAS    4.84    CM. 
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